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ABSTRACT 


This  thesis  describes  various  aspects  of 
the  origin  and  nature  of  microseisms,  in  the  light  of 
experimental  work  and  measurements  made  in  Alberta  during 
1962  and  1963. 

Microseisms  of  two  different  ranges  are  studied: 
3  to  10  seconds  and  2  to  50  cycles  per  second.  The  latter 
range  is  divided  into  the  3  cps  anomalous  peak  and  the 
4-50  cps  microseisms. 

The  3-10  second  microseisms  were  recorded  at 
the  University  of  Alberta  Seismic  Observatory  at  Edmonton 
during  July  and  August  1963-  Directions  of  approach  are 
determined  for  microseismic  storms  and  they  are  compared 
with  the  weather  situation.  It  is  concluded  that  3-5 
second  microseisms  of  Pacific  origin  are  very  heavily 
damped,  so  that  a  barrier  for  such  microseisms  exists. 

The  barrier  will  possibly  also  act  on  5  -  6  second  micro- 
seisms,  but  the  8  second  microseisms  of  Pacific  origin 
are  well  recorded. 

The  general  agreement  between  the  periods 
observed  and  those  predicted  from  Longuet-Higgins 1  theory 
adds  more  support  to  this  theory,  which  is  accepted  by 
some  but  not  by  all  seismologists. 
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A  peak  of  energy  at  3  cps  is  often  found  in 
Alberta,  and  it  is  shown  that  trains,  heavy  trucks,  and 
farm  machinery  are  effective  sources  for  such  microseisms. 

The  attenuation  and  the  nature  of  the  waves  are 
studied.  In  some  cases  retrograde  motion  is  dominant; 
in  other  cases  pure  prograde  motion  occurs  consistently. 
This  suggests  the.  occurrence  of  both  M]_  and  M2  waves  in 
the  3  cps  microseisms. 

A  clear  daily  variation  of  2-3  cps  microseisms 
is  demonstrated  and  it  is  shown,  that  this  is  not  due  to 
variation  of  wind;  other  causes  are  discussed. 

Microseisms  of  frequencies  7  and  12  cps  show 
consistent  prograde  particle  motion  in  some  cases,  but 
consistent  retrograde  particle  motion  in  others. 

In  the  Appendix  a  general  discussion  of  Rayleigh 
type  waves  in  a  two  layer  case  is  given.  It  is  shown  that 
Rayleigh's  constant  can  be  used  as  an  exploration  tool, 
as  suggested  by  Lee. 
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INTRODUCTION 


Historical,  Theories  on  Micro seisms 


Microseisms  are  defined  as  more  or  less  regular 
ground  motion  with  periods  not  exceeding  a  few  minutes  and 
continuing  with  constant  or  varying  amplitude  in  such  a  way 
that  they  can  be  assumed  not  to  be  caused  by  earthquakes  or 
explosions.  Exploration  geophysicists  use  the  term  ’’seismic 
noise"  to  designate  the  high  frequency  disturbances  encountered 
in  their  work.  They,  however,  include  the  shot- generated 
surface  waves,  when  they  are  unwanted  disturbances  on  the 
seismogram.  Since  the  present  work  deals  with  ground  vibra¬ 
tions  generated  by  sources  other  than  earthquakes  and  explos¬ 
ions,  the  term  "microseisms"  is  preferred. 

In  the  nineteenth  century  microseisms  were  noticed 
by  astronomers  and  geodesists,  who,  using  a  pool  of  mercury, 
rarely  found  the  surface  quiet.  And  Darwin  trying  to  observe 
the  lunar  tides  in  the  earth’s  crust  had  to  abandon  the 
experiment  because  of  ground  vibrations. 

The  term  "microseism"  was  introduced  by  Bertelli 
about  1870 ;  he  was  the  first  to  study  the  phenomenon  system¬ 
atically,  measuring  pendulum  movements  by  microscope;  and  he 
found  that  a  relationship  existed  between  microseisms  and 
disturbed  air-pressure  at  sea. 

After  the  improvements  of  seismic  instruments 
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beginning  late  in  the  nineteenth  century,  microseisms  attracted 
the  interest  of  many  observers.  Wiechert  (1905)  suggested 
the  surf  to  be  the  cause  of  the  observed  microseisms,,  and  was 
supported  by  Gutenberg  (1912  and  1921)  who  found  a  strong 
correlation  between  microseisms  in  Gottingen,  and  waves  observed 
at  Norwegian  coast  stations.  He  also  found  that  these  .micro- 
seisms  may  be  traced  as  far  away  as  Siberia. 

Through  many  years  the  surf  theory  was  the  only 
accepted  theory  for  generation  of  3-10  second  microseisms; 
in  German  it  was  even  called  "Brandungsunruhe " .  However 
Banerji  (1930)  noted  that  in  the  case  of  a  storm  approaching 
the  Indian  Ocean,  the  largest  microseisms  occurred  some 
hours  before  the  storm  reached  the  coastline.  The  same  result 
was  obtained  by  Gherzi  (1924),  who  put  forward  the  theory 
that  microseisms  are  generated  by  rapid  pressure  variations 
in  the  center  of  a  cyclone.  Several  objections  have  been 
made  to  this  theory:  the  amplitude  of  the  barometric  oscil¬ 
lations  seems  too  small,  and  the  theory  does  not  explain 
why  the  cyclone  has  to  be  situated  over  the  sea,  and  why 
the  sea  has  to  be  open  (not  covered  with  ice) . 

The  latter  conditions  are  satisfied  by  the  surf 
theory,  and  to  meet  Banerji1 s  objection  Gutenberg  and  van 
Straten  (1953)  use  the  term  "surf  theory"  in  a  broader  sense. 
This  means  that  the  new  surf  theory  only  states  that  the 
wave  energy  produced  by  storms  radiates  through  the  ocean 
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in  the  form  of  ocean  swell  until  it  reaches  a  continental 
boundary  or  the  coast.,  and  here  by  a  mechanism  not  specified 
in  the  theory,,  the  disturbance  is  transmitted  to  the  solid 
continental  block. 

Press  and  Ewing  (1948)  worked  out  a  mathematical 
theory  for  the  mechanism  of  the  energy  transmission  ard. 
found  that  the  periods  observed  in  the  ordinary  microseisms 
(i.e.  3-10  sec.)  could  be  explained,  regarding  the  system 
water  plus  bottom  as  a  single  acoustic  system  giving  normal 
mode  propagation  with  constructive  interference,  when  periods 
and  water  depths  were  subject  to  certain  conditions. 

How  the  oscillations  are  excited  is  not  quantita¬ 
tively  specified  in  Press  and  Ewing’s  theory.  Donn  (1957) 
proposes  that  the  breaking  of  waves  at  sea  would  produce  a 
sufficiently  high  random  noise  level,  while  Longuet-Higgins 
(1950)  develops  a  mathematical  theory,  based  on  Miche’s 
finding  that  second  order  pressure  variations  due  to  inter¬ 
ference  of  opposing  waves  can  produce  microseisms  of  the 
required  amplitude.  The  theory  predicts  a  2:1  ratio  of  water 
wave  periods  and  microseism  periods  and  often  such  ratios 
have  been  found.  The  theory  has  subsequently  been  verified 
in  tank  experiments.  Banerji  (1935)  shows  by  tank  experi¬ 
ments  that  pressure  variations  can  be  propagated  with  acoustic 
speed  from  a  certain  plane  near  the  surface  =a=ndr  to  the  bottom. 
In  most  treatises  on  water  waves  these  pressure  variations  are 
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ignored,  because  they  do  not  appear  when  the  usual  procedure 
of  dropping  second  order  terms  is  applied. 

The  assertion  that  almost  all  microseisms  in  the 
period  spectrum  3-8  seconds  are  generated  by  the  mechanism 
of  interference  of  water  waves  is  believed  by  some  and 
questioned  by  others,  but  it  may  be  regarded  a  good  working 
hypothesis . 

Haubrich  et  al  (1963):,  using  a  digital  seismograph 
system,  find  a  seismic  spectral  peak  with  a  frequency  which 
is  twice  the  frequency  of  the  ocean  waves  as  3=^  is  expected 
according  to  Longuet-Higgins 1  theory,  but  they  also  find  a 
seismic  peak  of  the  same  frequency  as  the  ocean  wave  frequency, 
as  expected  according  to  Wiechert's  theory  (the  frequency 
quoted  is  70  millicycles  per  second) .  The  latter  seismic 
peak  has  only  one  per  cent  of  the  power  which  is  associated 
with  the  first  one. 

There  are  different  reasons  why  the  microseisms 
are  studied  so  eagerly;  one  is  that  microseisms  show  many 
unexplained  features  like  narrow  spectra,  regular  beats  or 
rapidly  varying  amplitudes.  Another  reason  is  that  the  study 
of  microseisms  might  make  it  possible  to  infer  something 
about  parameters  of  the  earth's  crust  and  about  atmospheric 
and  oceanographic  disturbances.  Very  often  conclusions  of 
this  kind  have  been  made  in  publications  concerning  a  few 
case  histories  but  very  seldom  have  the  results  achieved 
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been  generally  accepted. 

The  above  discussion  Is  concerned  with  the  3-10 
second  microseisms.  In  recent  years,,  however,  the  interest 
in  microseisms  of  higher  frequencies  has  been  increasing. 
This  is  partly  because  of  the  more  widespread  use  of  instru¬ 
ments  capable  of  recording  higher  frequencies,  and  partly 
because  of  the  importance  of  choosing  quiet  sites  for  the 
costly,  but  efficient  observatories  using  seismic  arrays. 
Moreover,  "the  study  of  short  period  microseisms  could 
become  of  interest  to  the  earthquake  engineers,  because 
E.  M.  Antonenko  (1963)  in  her  studies  of  microseisms  found 
a  close  correlation  between  the  amount  of  old  earthquake 
damage  and  the  amplitude  of  microseisms  in  the  period  range 
0.15  to  0.35  seconds",  (Y.  F.  Savarensky,  personal  communi¬ 
cation)  . 

During  the  last  fifteen  years,  only  13  out  of  238 
papers  on  microseisms,  abstracted  by  the  U.  S.  Geological 
Survey  Bulletins^  were  concerned  with  short  period  (less 
than  one  second)  microseisms,  but  in  future  this  ratio  is 
likely  to  increase.  Suggested  sources  for  these  microseisms 
are  traffic,  machinery,  wind  and  rain  (Wilson,  1953a) >  atmo¬ 
spheric  micro-oscillations  (Walsh,  1955)  and  daylight 
(Gherzi,  1962).  Further  discussions  on  short  period  micro¬ 
seisms  follow  in  chapter  5,  together  with  some  results 
based  on  the  University  of  Alberta  recordings. 
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The  main  object  of  the  approach  made  in  this  thesis 
is  to  investigate  microseism  sources  and  the  nature  of  the 
waves.  Another  approach  is  made  by  Galbraith  (1963)  who 
treats  microseisms  from  a  statistical  point  of  view  for  the 
purpose  of  designing  filters  that  would  suppress  the  micro¬ 
seisms  and  enhance  the  seismic  signals.  The  approach  in  this 
thesis  is  to  do  something  with  them;  Galbraith’s  approach  is 
to  do  something  about  them. 

1 . 2  Nature  of  Microseismic  Waves 

Numerous  authors  have  found  that  microseisms  contain 
a  large  fraction  of  ordinary  Rayleigh  waves*;  some  have  even 
come  to  the  conclusion  that  all  of  them  are  ordinary  Rayleigh 
waves,  while  in  other  cases  comparable  proportions  of  Rayleigh 

SL  1L. 

and  Love  waves  '  have  been  demonstrated.  Thus  it  appears  that 
the  vertical  component  of  microseisms  has  often  been  assumed 
to  consist  entirely  of  Rayleigh  waves  and  it  is  only  close 
to  the  source  (Wilson  1953b)  that  body  waves  have  been 
demonstrated  as  an  additional  component.  This  of  course  does 
not  exclude  the  possibility  that  other  kinds  of  surface  waves 
than  ordinary  Rayleigh  and  Love  waves  could  be  components  of 
microseisms . 

One  example  demonstrating  this  possibility  was 
given  by  Monakhov  and  Dolbilkina  (1958)  in  an  investigation 
of  microseisms  at  Yalta.  Using  an  ordinary  seismograph,  an 

Surface  waves  with  a  longitudinal  and  a  vertical  transverse 
displacement . 

**  Surface  waves  with  a  horizontal  transverse  displacement. 

***  Elastic  waves,  either  longitudinal  P  waves  or  transverse 
S  waves . 
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azimuthal  seismograph,,  a  vector  seismograph  and  a  tripartite 
station.  When  studying  individual  oscillations  of  Atlantic 
origin,,  they  found  a  complicated  particle  motion  in  most  of 
the  cases,  hut  a  narrow  angle  of  approach  was  usually 
pointing  towards  the  Norwegian  coast.  In  seven  to  twenty- 
seven  per  cent  of  the  cases  they  found  pseudo-Rayleigh  waves, 
i.e.,  waves  differing  from  pure  Rayleigh  waves  only  by  having 
a  plane  of  oscillations  tilted  towards  one  side,  but  retaining 
all  other  properties  of  the  Rayleigh  waves;  they  found  prac¬ 
tically  no  Love  waves.  It  is  tempting  to  suggest  that  those 
pseudo-Rayleigh  waves  are  due  to  interference  between  pure 
Rayleigh  waves  and  Love  waves  because  when  only  a  few  consecu¬ 
tive  oscillations  can  be  studied  it  seems  impossible  to 
distinguish  between  the  two  phenomena.  However,  if  Love 
waves  were  present  they  should  sometimes  occur  between  the 
beats  of  the  vertical,  and  the  azimuthal  setup  would  make 
their  identification  probable;  but  no  such  Love  waves  were 
reported.  Furthermore,  in  a  large  number  of  samples,  the 
Love  waves  would  be  statistically  independent  of  the  Rayleigh 
waves  and  the  apparent  angle  of  inclination  should  occur  on 
both  sides  of  the  vertical  with  equal  probability.  However, 
for  most  azimuths  there  appeared  to  be  a  strong  tendency  for 
the  plane  of  vibration  always  to  tilt  to  the  same  side.  The 
angle  of  tilt  was  most  often  twenty  to  thirty  degrees  and 
never  exceeded  forty-five  degrees. 

This  phenomenon  seems  unexplainable,  using  simple 
horizontal  layer  models,  but  may  be  explainable  using  a  more 
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complicated  crustal  model  for  southern  Russia. 

An  interpretation  of  these  data  without  recourse 
to  departures  from  classical  elasticity  theory  can  be  derived 
from  the  discussions  given  by  Oliver  and  Ewing  (1958a)  and 
Jones  (1963):  They  describe  Leet’s  C  wave  or  coupled  wave, 
which  has  a  definition  similar  to  that  of  the  pseudo-Rayleigh 
wave.  They  point  out  that  Love  waves  and  Rayleigh  waves  would 
arrive  almost  simultaneously  if  the  travel  times  of  such 
waves  are  similar,  and  they  show  that  a  similarity  exists 
between  the  velocities  of  continental  Rayleigh  waves  and 
of  the  second  mode  Love  wave. 

Since  pseudo-Rayleigh  waves  were  found  to  have  a 
fairly  small  angle  of  inclination  at  Yalta,  the  transverse 
component  of  such  waves  is  expected  to  be  smaller  than  the 
longitudinal.  However,  Jensen  (1961)  found  no  significant 
difference  between  the  amplitudes  cf  the  longitudinal  and 
transverse  component  for  those  eighty-two  days  in  which 
the  microseisms  were  approaching  from  the  northern  octant. 

This  strongly  suggests  that  Love  waves  are  present  in  the 
microseisms  recorded  at  Copenhagen,  and  that  their  amplitudes 
are  similar  to  those  of  the  Rayleigh  waves.  Furthermore,  the 
sinusoidal  variation  of  velocity  vectors  (Jensen  1958,  p.  16) 
are  most  likely  due  to  interference  of  a  group  of  Love  waves 
and  a  group  of  Rayleigh  waves.  It  seems  improbable  that 
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two  simultaneous  but  entirely  different  sources  should  make 
a  regular  interference  pattern  such  as  was  suggested  by 
Bath  (1962). 

The  Rayleigh  waves  and  the  pseudo-Rayleigh  waves 
discussed  above  both  have  a  retrograde  particle  motion, 
i.e.  motion  with  a  velocity  towards  the  source  in  the  upper 
portion  of  the  trajectory.  The  opposite  type  of  particle 
motion  is  called  prograde,  and  this  particle  motion  is  the 
same  as  that  exhibited  by  water  waves. 

Prograde  waves  were  theoretically  predicted  by 
Sezawa  and  Kanai  (1935)  by  considering  the  dispersion  of 
Rayleigh  waves  in  a  two  layer  case.  The  equations  resulted 
in  two  different  dispersion  curves,  one  of  them  corresponded 
to  the  fundamental  Rayleigh  mode,  the  other  one  they  called 
the  Mg-mode,  both  modes  were  considered  fundamental  because 
of  the  absence  of  nodal  planes.  They  found  that  the  M2- 
Rayleigh  wave  (Sezawa  wave)  had  orbital  motion  opposite 
that  of  ordinary  Rayleigh  waves  (i.e.,  prograde  instead  of 
retrograde).  Nagamune  (1956)  found  that  in  a  three  layer 
case  M2  waves  may  be  either  prograde  or  retrograde. 

Oliver  and  Ewing  (1957)  pointed  out  that  the 
calculations  of  Sezawa  and  Kanai  had  led  to  the  general 
belief  that  the  M2  mode  was  characterized  by  prograde 
(elliptical  and  progressive)  particle  motion.  In  a  two 
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layer  case  with  the  low  rigidity  contrast  (1:2.09)  pertinent 
to  the  problem  of  higher  mode  continental  Rayleigh  waves, 
however,  they  found  the  wave  to  have  retrograde  motion  for 
all  wave-lengths,  and  in  another  case  with  higher  rigidity 
contrast  (1:8),  they  found  the  particle  motion  to  change 
from  prograde  to  retrograde  within  the  same  mode . 

While  M2  waves  may  have  both  types  of  motion,  the 
fundamental  mode  Rayleigh  wave,  at  least  in  any  realistic, 
layered  model,  will  always  have  retrograde  particle  motion 
(J.  Oliver,  personal  communication). 

In  1955  Kanai  ascertained  that  M2  waves  appear 
only  when  the  velocity  in  the  surface  layer  is  smaller  or 
slightly  larger  than  the  velocity  in  the  lower  layer  and  he 
concluded : 

"Moreover,  though  Sezawa  waves  cannot  be  seen 
when  the  wave-length  exceeds  a  certain  limit,  the  reason 
why  Sezawa  waves  are  rarer  than  Rayleigh  waves  in  cases  of 
shorter  wave-lengths  is  hard  to  find.  Also,  there  is  no 
reason  why  the  amplitude  of  Sezawa  waves  is  smaller  than 
that  of  Rayleigh  waves,  but  in  some  cases,  it  seems  that 
the  amplitude  of  Sezawa  waves  becomes  larger  than  that  of 
Rayleigh  type  waves." 

In  view  of  this  theoretical  prediction  it  is  interesting  to 
note  that  prograde  Rayleigh  waves  later  were  found  in  waves 
from  explosions,  as  demonstrated  by  Kisslinger  (1959)  and 
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by  Jones  (1963).  They  have  also  been  found  In  microseisms,, 
which  will  be  demonstraded  in  the  present  work.  No  previous 
demonstration  of  such  Sezawa-type  waves  in  microseisms  is 
known  to  the  writer.  It  is  not  very  surprising,  though,  that 
previous  workers  have  felt  justified  in  the  assumption  that 
the  vertical  component  of  microseisms  was  composed  of  retro¬ 
grade  Rayleigh  waves.  The  reason  is  that  earlier  work  on 
particle  motion  in  microseisms  has  been  done  on  microseisms 
with  periods  longer  than  four  seconds,  and  several  investi¬ 
gations  suggest*  that  any  prograde  waves  in  such  microseisms 
are  negligible.  One  such  investigation  is  the  routine  deter¬ 
mination  of  the  direction  of  approach  in  Copenhagen  by  Jensen’s 
method  (Jensen  1959).  Here  the  most  frequent  direction  of 
approach  was  from  the  northern  sector,  confirming  Gutenberg’s 
original  work  and  recent  Russian  work  indicating  the  Norwegian 
coast  as  a  prominent  source  area.  However,  if  the  microseisms 
had  contained  an  appreciable  amount  of  prograde  waves,  a 
number  of  directions  should  group  around  the  antipole,  i.e., 
in  the  southern  sector,  and  no  such  tendency  occurred. 

Another  investigation  made  by  Oliver  and  Ewing 
(1957)  led  to  the  same  result.  They  demonstrated  theoretically 
(as  mentioned  before)  and  experimentally  (investigating  an 
Arctic  earthquake)  that  the  M2  particle  motion  is  retro¬ 
grade  in  a  continental  crust.  The  microseisms  in  which 
prograde  particle  motions  were  observed  at  the 
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University  of  Alberta  were  of  periods  shorter  than  half  a 
second . 

Some  additional  Information  about  Rayleigh- type 
waves  In  a  two  layer  model  Is  given  In  Appendix  A. 
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II.  METHODS  OF  INVESTIGATIONS 

2 . 1  Procedure  for  Investigating  the  Nature  of  Microseismic 

Waves 

Jensen’s  method  (see  page  28)  is  a  method  for  deter¬ 
mining  the  direction  of  approach,  given  the  assumption  that 
Rayleigh  waves  are  present;  the  same  method  will  be  used  here 
for  determining  the  nature  of  waves  assuming  that  the  direction 
of  approach  is  known.  Since  the  purpose  thus  is  modified,  it 
could  be  called  the  modified  Jensen  method. 

The  first  concern  of  this  method  is  knowledge  of 
the  source.  In  some  cases,  particularly  in  artificial  noise, 
the  source  can  be  established  beyond  any  doubt.  In  such 
cases  the  direction  of  approach  will  most  often  be  the  direc¬ 
tion  to  the  source,  but  some  lateral  refraction  may  occur 
on  the  path.  A  safer  method  is  therefore  to  use  an  L-spread 
of  geophones  to  determine  the  direction  of  approach  right  on 
the  site,  where  the  waves  are  analyzed.  In  addition  to  this 
a  three-dimensional  registration  must  be  made:  one  vertical 
(Z),  and  two  horizontals,  usually  N  and  E,  all  recorded  at 
the  same  location.  The  three  components  should  preferably 
have  the  same  responses.  If  that  is  the  case  the  procedure 
is  to  select  times  where  Z  is  maximum  and  measure  dN/dt  and 
dE/dt  and  plot  these  points  in  a  NE  diagram.  In  the  same 
diagram  the  vector  pointing  towards  the  source  and  the  line 
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through  the  origin  perpendicular  to  this  line  are  drawn. 

Now  the  distribution  of  points  obtained  is  compared 
with  the  cases  shown  in  Figure  1,  and  it  may  be  possible  to 
decide  which  case  is  appropriate.  If  small  amounts  of 
extraneous  microseisms  are  present  the  scatter  may  be  increased 
somewhat,  and  if  larger  amounts  are  present  it  may  be  necessary 
to  use  statistical  methods  on  a  large  number  of  determinations 
before  any  conclusions  can  be  drawn,  as  shown  by  Jensen  (1961). 
If  in  cases  (b)  and  (c)  Rayleigh  waves  and  Love  waves  have 
constant  amplitudes  of  similar  magnitudes  the  points  would 
scatter  in  a  quarter  plane,  but  due  to  the  amplitude  modula¬ 
tions  usually  present  in  microseisms,  they  will  scatter  in 
a  half  plane  as  shown  in  Figure  1. 

2 . 2  Procedure  for  Determination  of  the  Direction  of  Approach 
by  Wave  Correlation  in  an  L-Spread 

When  three  or  more  suitably  spaced  geophones  are 
laid  out  to  form  a  right  angle  or  a  cross,  it  is  possible 
to  determine  the  direction  of  approach,  provided  that  a  sub¬ 
stantial  part  of  the  waves  in  a  fairly  narrow  frequency 
range  have  their  directions  of  approach  in  a  fairly  narrow 
sector.  Successes  of  the  tripartite  method  have  been  obtained 
when  this  condition  was  satisfied,  and  failures  have  resulted 
in  cases  when  it  was  not  satisfied.  When  the  condition  is 
not  satisfied,  it  is  sometimes  possible  to  identify  waves 
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Possible  nature  of  surface  waves,  when  only  one  source  is  present: 

(a)  Mj  .Rayleigh  waves  (  ordinary  Rayleigh  waves) 

(b)  M|  Rayleigh  waves  and  Love  waves 

(c)  M2  RayJeigh  waves  and  Love  waves 

(d)  M2  Rayleigh  wa  ves  (  Sez  awa  waves  ) 

Points  plotted  are(^' rea^  Qt  times  where  Z  is  maximum 

(  artificial  example  ) 
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from  different  directions,,  and  this  can  he  done  with  much 
better  chances  of  success  when  a  spread  of  many  geophones 
is  used  in  each  direction  instead  of  just  two  as  used  in 
the  tripartite  method.  The  advantage  of  having  more  than 
two  geophones  is  that  a  propagating  wave  crest  or  trough 
will  show  a  clear  lineup,  which  often  can  be  identified  on 
top  of  other  disturbances.  One  example  of  a  mathematical 
evaluation  of  the  tripartite  method  was  given  by  Gutenberg 
(1958)  for  an  arbitrary  location  of  the  three  stations.  In 
the  text  below  a  simpler  method  is  evaluated,  which  applies 
when  it  is  possible  to  determine  the  apparent  velocity  in  two 
mutually  perpendicular  directions. 

Assume  apparent  velocities  and  along  a  NS 
line  and  an  EW  line  are  known.  Positive  directions  are  taken 

to  be  N  and  E.  Consider  a  wave  propagating  in  a  direction 

A  I  A  I  A  A 

determined  by  the  vector  v.  Let  I v I  =  v.  Let  n  and  e  be 

unit  vectors  pointing  toward  N  and  E.  The  angle  between  n 

and  (-v)  is  called  a,  i.e.,  a  is  the  azimuth  of  the  direction 

of  approach.  Let  a  be  positive  in  the  two  eastern  quadrants 

and  negative  in  the  two  western  quadrants,  -l80°  <  a  <:  l80° . 

A 

Consider  a  point  moving  with  the  wave  in  the  direction  v. 

A  A 

The  velocity  of  its  projection  on  a  meridian  is  then  v  •  n, 
positive  if  moving  north,  and  the  velocity  of  its  projection 
on  a  parallel  is  v  •  e  positive  if  moving  east.  v^  has  the 
same  sign  as  v  •  n  and  hence  (Figure  2) 
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ol  :  Direction  of  approach 
v  :  Direction  of  wave  .propagation 
v=lvl  :  Velocity  of  wave  propagation 
vN  :  Apparent  velocity  in  SN-direciion 
vfc-  :  Apparent  velocity  in  W£~direction 

Determination  of  the  direction  of  approach  from  the  apparent 

wave  velocity  in  an  L  spread. 


FIG.  2. 
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VE  has  the  same  sign  as  v  •  e  and  hence  (Figure  2) 
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from  which  it  follows  that 
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As  shown  by  Longuet-Higgins  (1957)  the  relations  are 


simpler  when  vectorial  wave  numbers  are  used. 
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«  =  WM  is  the  wave  number  in  the  N-direction 
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Wave  number  in  a  direction  given  by  the  unit  vector  m  is  w  •  m 
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If  several  waves  are  present  simultaneously  with 
wave  numbers  0Pn  which  are  slightly  different,,  then  the  wave 
number  determined  from  the  L-spread  will  be  the  center  of  mass 
of  the  points  Pn,  each  weighted  by  the  corresponding  component 
amplitude.  In  that  case  the  velocity  determined  is  the  aver¬ 
age  phase  velocity,  but  as  was  pointed  out  by  Jones  (1963 
and  personal  communication)  the  differences  between  the  com¬ 
ponent  wave  numbers  (and  hence  frequencies)  must  be  so  small 
that  their  differences  only  contribute  to  the  shape  of  the 
envelope  (the  modulation)  without  too  much  distortion  of  the 
carrier  frequency. 

2 . 3  Procedure  for  Measurement  of  Amplitudes  and  Frequencies 

of  Microseisms 

Most  often  microseisms  show  a  beat  phenomenon,  and 
irregular  microseisms  are  often  made  up  of  different  more 
regular  components  each  of  which  shows  beats,  so  that  each 
component  frequency  may  show  up  as  a  sinusoidal  train  of  waves 
between  beats  of  other  frequencies.  It  is  thus  possible  to 
measure  some  component  frequencies  by  visual  observation,  and 
if  the  recording  is  not  made  on  magnetic  tape,  suitable  for 
automatic  power  spectrum  analysis,  a  fairly  accurate  measure 
of  the  power  associated  with  a  certain  frequency  can  be 
obtained  by  taking  the  average  of  several  typical  beat  maxi¬ 
mum  amplitudes.  This  yields  an  estimate  somewhat  larger  than 


the  average  amplitude  (in  the  mean  square  sense),  but  it  is 
an  estimate  that  according  to  practical  experience  and 
theoretical  considerations  doesn’t  vary  too  much  as  long  as 
the  source  conditions  remain  constant. 

Such  theoretical  considerations  were  given  by  Longuet- 
Higgins  (Vesiac  Staff,  1962,  Chapter  8) ,  who  defines  a  para¬ 
meter  e  related  to  the  narrowness  of  the  spectrum.  An  esti¬ 
mate  of  €  is  given  by 


e 


2 


(2-1) 


where  N  is  the  number  of  zero  up  crossings 
o 

Nq1  is  the  number  of  crests 
in  some  length  of  the  record. 

When  €  =  o  the  median  of  the  distribution  of  the 
maximum  crest  amplitudes  in  a  random  sample  of  N  crests  will 
be  2.9  times  the  root  mean  square  amplitude  at  N  =  64  and 
3.4  times  rms  amplitude  at  N  =  256.  The  corresponding  10 
percent  fractiles  are  2.6  and  3.1,  i.e.,  the  ratio  between 
maximum  crest  and  rms  amplitude  will  usually  be  around  3. 

This  is  for  gaussian  noise  with  e  =  o,  but  also  for  gaussian 
noise  with  e  =  0.9  the  ratio  will  be  around  3  for  a  reasonably- 
sized  sample.  A  reason  for  not  choosing  larger  samples  is 
that  the  condition  of  statistical  stationarity  is  not  often 
fulfilled  in  a  sample  larger  than  a  few  hundred. 
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Figure  3  shows  the  cumulative  probability  distri¬ 
butions  of  a  quantity  x,  which  is  the  maximum  crest  amplitude 
divided  by  the  rms  amplitude,  in  a  sample  of  N  crests ,  for 

N  =  2°,  2^,  . ,  2 (taken  from  Cartwright,  1958). 

The  value  obtained  by  this  procedure  will  thus  be 
approximately  three  times  the  rms  amplitude  for  gaussian 
noise.  When  gaussian  noise  has  a  prominent  and  fairly  narrow 
spectrum  (most  of  the  energy  within  one  octave)  it  is  possible 
by  visual  observation  to  obtain  the  mean  frequency  of  this 
power  density  peak.  This  was  shown  by  Rice  (1944,  1954)  who 
proved  that  the  mean  number  of  zeros  per  second  is 

2  [  ^  ]l/2  (2-2) 

mQ 

where  the  moments mR  of  the  spectrum  are  defined: 

mn  =  /“  w(f )fndf  (2-3) 

w(f)  being  the  power  density  associated  with  the  frequency  f. 
For  n  =  0 

m0  =  r  w(f)df  =  11m  fl  [y(t)]2  dt  (2-4) 

o  o  T_oo  21  -1 

where  y(t)  is  the  ground  amplitude  measured  from  its  mean 
level.  /mQ  is  the  rms  value  of  the  amplitudes. 

The  interpretation  of  w(f)  as  power  density  requires 
a  little  clarification.  For  frequencies  between  f  and  f  +  df 
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Figures  1  a  and  b.  Probability  functions  Pn{%)  (equations  (1)  and  (4))  for  various  values  of 

N  =  2m  and  (a)  e  =  0,  (b)  e  =  0-9. 

After  Cartwright  ( 1958 ) 

FIG. 3. 

Probability  functions  for  the  amplitude  of  the  highest 


of  N  crests. 
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the  power  output  of  the  seismometer  will  he  proportional 
to  w(f)df,  with  the  same  constant  of  proportionality  for 
any  location.  Thus  w(f)  is  a  measure  of  the  power  density 
of  seismometer  output.  For  stationary  surface  waves  w(f)df 
and  hence  mQ  will  he  proportional  to  the  mean  energy  per 
unit  area  of  the  surface,,  hut  the  constant  of  proportionality 
will  depend  on  the  properties  of  the  ground.  Thus  a  larger 
power  density  on  one  location  than  another  does  not  neces¬ 
sarily  mean  larger  ground  energy;  if  the  ground  is  softer  it 
may  even  he  lower. 

Longuet-Higgins  (Vesiac  staff,  1962,  Chapter  8) 
pointed  out  that  in  gaussian  noise  the  distrihution  of 
intervals  t  between  successive  zeros  has  the  following  proper¬ 
ties  when  the  narrow  hand  approximation  holds: 

a)  It  is  symmetrical  about  r  =  -i(— 

d  m.r) 

h)  It  is  independent  of  the  shape  of  the  original 
spectrum  except  for  its  rms  width  6 
The  narrow  hand  approximation  is  valid  when  6  «  1 

where : 

&2  _  m2m0  ~  ml  _  (f2)  -  (f)2  =  (f  -  f)2 

m2m0  Iff  Iff 

(2-5) 

it  holds  when  most  of  the  energy  is  within  one  octave  or 


less . 


? 
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If  the  energy  is  divided  between  two  or  more  power 
density  peaks  more  than  an  octave  apart,,  the  shorter  periods 
will  clearly  show  up  superimposed  on  the  longer  ones,  and  each 
kind  of  microseism  may  be  analyzed  separately  from  the  same 
record . 

Thus  for  an  analogue  record  it  is  possible  to  deter¬ 
mine  the  mean  frequency  by  simply  taking  the  mean  of  zero- 
crossover  intervals  at  randomly-distributed  times.  When  e  =  0, 
which  is  often  the  case,  the  interval  between  successive  crests 
will  simply  be  twice  the  interval  between  successive  zeros. 

When  €  >  0,  the  distribution  of  crests  will  involve  the  moment 
of  fourth  order  (Vesiac  Staff,  1962,  Chapter  8) .  According 
to  the  same  source  the  empirical  distribution  of  zero- crossover 
intervals  will  also  depend  on  the  width  of  the  spectrum,  so 
that  the  interquartile  range  of  this  distribution  will  be 
equal  to  [2//3](  6  t)  .  By  estimating  the  interquartile  range 
it  is  thus  possible  to  form  an  estimate  of  6,  the  rms  width 
of  the  spectrum  per  unit  of  frequency. 

If  more  knowledge  of  the  shape  of  a  spectral  peak 
is  required,  or  if  resolution  of  two  closely-spaced  spectral 
peaks  is  essential,  digitizing  the  record  is  required,  in 
order  to  calculate  the  power  spectrum. 


’ 

* 


25. 


2.4  Procedure  for  Determination  of  the  Attenuation  of  Micro¬ 


seisms 


When  the  term  "attenuation  of  seismic  waves"  is 


used,  it  means  the  decrease  of  amplitude,  which  theoretically 
is  due  to  conversion  of  elastic  energy  to  heat,  but  in  prac¬ 
tice  scattering  and  reflections  may  account  for  a  portion  of 
the  attenuation.  At  any  rate  the  attenuation  will  be  zero, 
in  a  perfectly  elastic,  homogeneous  medium.  Bullen  (1953) 
prefers  the  word  damping,  meaning  precisely  the  same  as  what 
other  authors  would  call  attenuation. 

Since  microseisms  consist  primarily  of  surface 
waves,  it  is  reasonable  to  define  attenuation  of  microseisms 


with  reference  to  Rayleigh  waves  in  a  homogeneous,  isotropic, 
perfectly  elastic  half space .  For  such  a  half space  Rayleigh 
waves  from  an  infinite  line  source  would  have  a  constant 
amplitude  as  a  function  of  distance.  Rayleigh  waves  from  a 


io)t 


point  source,  with  harmonic  motion  proportional  to  e 


however,  would  have  the  following  vertical  amplitude  (Ewing, 
Jardetzky  and  Press,  1957  j  page  40)  : 


w 


o 


(2-6) 


where 
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and  the  constants  a,  f3,  CU  are  the  velocities  of  P,  S  and 

n 

Rayleigh  waves  respectively,  and  a)  =  2irf . 

The  Bessel  function  in  (2-6)  is  the  only  factor 
that  depends  on  distance.  When  kr  is  large  enough  JQ(kr) 
will  approximate  a  sine  wave  with  constant  wave  number  k. 
This  result  follows  from  the  approximation  formula  (Gray  and 
Mathews,  1895) : 

J0U)  -  (^)1//2  M  cos  (x  -  f  ■  (2-7) 

where 

1  58 

M  =  1  -  — "  . 

l6x2  512x4 

ip  =  tan'1  (-gj  -  ^21^  +  .  .  .  .) 

It  also  follows  that  when  r  is  large  the  amplitude 
of  the  Rayleigh  wave  will  be: 

A  =  Br"1/2  (2-8) 

where  the  factor  B  =  Ck^2  vR  • is  independent  of  r. 

In  a  medium  where  attentuation  occurs  (2-8)  will 
be  modified  into: 

A  =  B  •  r'1/2  •  e“ar  '  (2-9) 

and  a  may  be  called  the  attenuation  coefficient.  Kisslinger 
(Jones,  1963)  found  a  to  vary  from  4  to  6  km”1  in  the  range 
3.5  to  5  cps.  Wilson  (1953b)  found  a  =  1.09  km”1  at  9  cps . 
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Knopoff  and  McDonald  (1958)  defined  the  specific 
dissipation  function  l/Q  such  that  2tt/Q  is  the  ratio  of 
energy  dissipated  per  cycle  to  the  energy  stored.  For  a 
plane  propagating  elastic  wave  they  found: 


where  c  is  the  velocity  of  propagation. 

Knopoff  (1959)  quoted  two  different  models  for  the 
attenuation : 

1.  The  "Q"  or  solid  friction  .model,  that  predicts  Q  to  be 
independent  of  the  frequency  f,  and  hence  the  attenuation 
coefficient  a  to  be  proportional  to  f. 

2.  The  internal  friction  model,  that  predicts  Q  to  be 

p 

proportional  to  f  and  a  to  be  proportional  to  f  . 

Knopoff  and  MacDonald  (1958)  quote  evidence  from 
many  sources  that  the  "Q"  model  applies  for  most  materials 
at  seismic  and  acoustic  frequencies,  and  it  may  be  expected 
to  apply  for  most  homogeneous  materials. 

Willis  and  Wilson  (i960)  found  =  140  -  152 
valid  for  1  cps  Rayleigh- type  waves  observed  100  -  1000  km 
from  explosive  sources.  This  agrees  with  the  range  Q  =  100  - 
250  given  by  Knopoff  (1959). 

Brune  (1962)  quoted  Q,  =  170  for  Rayleigh  waves 
with  period  T  =  140  -  215  seconds  and  Q  =  350  for 
T  =  250  -  350  seconds. 


He  concluded  that  for  dispersive 


’ 


28. 


wave  trains  the  equation  (2-9.1)  is  still  valid,  provided 
that  c  is  interpreted  as  the  group  velocity.  He  also  con¬ 
cluded  that  measured  values  of  Q  usually  reflect  the  combined 
loss  of  energy  due  to  internal  friction  and  scattering. 

If  an  attenuation  study  is  carried  out  using  a 
comparatively  small  number  of  observations,  it  is  necessary 
to  check  that  the  source  conditions  remain  constant.  Even  if 
the  source  is  an  explosion  of  the  same  size  in  the  same 
shot  hole,  the  energy  coupling  may  have  changed  at  the  source. 

Sometimes  the  surface  amplitude  of  seismic  signals 
shows  a  clear  variation  from  one  place  to  another.  For 
instance  earthquake  damage  will  be  larger  on  loose  sediments 
than  on  solid  rock.  Consequently  in  an  attenuation  study, 
it  must  be  considered  that  amplitude  differences  could  be 
due  to  such  an  effect. 

A  detailed  example  of  an  experimental  determination 
of  an  attenuation  coefficient  and  a  discussion  of  Q  values 
for  short  period  surface  waves  is  given  in  section  5 -3c. 

2 . 5  Jensen's  Method  for  the  Determination  of  the  Direction 

of  Approach  of  Rayleigh  Waves 

A  convenient  method  to  apply  on  photographic 
records  from  a  homogeneous  set  of  instruments  is  the  Jensen 
method.  The  first  part  of  this  section  will  contain  an 
outline  of  this  method,  while  the  later  parts  will  show  some 
new  theoretical  work,  which  the  author  found  pertinent  to  the 
interpretation  of  some  results  achieved  by  the  Jensen  method. 


■ 

■ 


29. 


2.5a  Short  Outline  of  Jensen* s  Method 

Jensen’s  method  or  the  method  of  the  empty  half 
plane  has  been  described  by  Jensen  (1958  and  1961)  using 
vectors  and  by  Bath  (1962)  using  Cartesian  coordinates. 

Basic  Assumptions: 

1)  Retrograde  Rayleigh  waves  from  a  single  direction  are 
present  and  they  dominate  over  any  other  wave type  with 
a  vertical  (Z)  component. 

2)  Any  other  wave types  contained  in  the  movement  are  statis¬ 
tically  independent  of  these  Rayleigh  waves  at  the  times 
when  Z  is  at  a  maximum. 

Theory: 

At  the  times  when  Z  is  at  a  maximum  the  retrograde  Rayleigh 
component  will  contribute  to  the  particle  velocity  by  a 
vector  vL  pointing  towards  the  source.  A  Love  wave  not 
correlated  with  the  Rayleigh  wave  will  contribute  to  the 
particle  velocity  by  a  vector  in  one  of  the  directions 
perpendicular  to  VR,  but  it  will  assume  these  two  directions 
with  equal  probability.  The  total  effect  of  these  two 
components  will  be  that  all  the  end  points  of  the  velocity 
vectors  scatter  in  the  half  plane  towards  the  source.  An 
addition  of  other  wavetypes  and  isotropic  noise,  uncorrelated 
with  the  Rayleigh  component  will  not  destroy  the  tendency  of 
these  points  to  choose  the  half  plane  towards  the  source. 
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Procedure : 

1.  Times  of  well-developed  maxima  on  Z  are  determined. 

dN  dE 

2.  The  slopes  and  are  measured  at  these  times. 

3.  The  points  (-^r ,  -^r)  are  plotted  in  an  NE-diagram. 

4.  Now  the  vector  (Z  -^r,  S  ■—)  is  an  estimate  of  the 
direction  of  approach. 

5.  The  perpendicular  to  this  vector  through  origin  is 
drawn  and  the  number  r  of  points  in  the  half  plane 
containing  the  vector  (the  "right"  half  plane)  and  the 
number  of  points  w  in  the  other  half  plane  (the  "wrong" 
half  plane)  are  counted. 

6.  The  reliability  measure  p  is  determined  by  the  formula 

p  =  •  10  (2-10) 

where  n  is  the  total  number  of  points.  An  approximate 
value  of  the  probable  error  of  the  direction  of  approach 
for  positive  p  values  is 

e  =  30  -  p  (2-11) 

where  e  is  the  probable  error  in  degrees  (Jensen  1961) . 

2.5b  Definition  and  Interpretation  of  Ap 

To  make  use  not  only  of  the  direction  of  the 

vector  (Z  ,  2  -^|)  quoted  in  the  previous  sub-section,  but 

also  of  its  absolute  value,  it  was  found  convenient  to  intro¬ 
duce  a  quantity  A^  defined  by 


, 
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Ts[(sH)2+  (S  ^|)2]1/2  •  0.798 


dx 


27m  M(T) 


where  T 


s 

x 

M(T) 


=  period 
=  paper  speed 

=  distance  along  recorded  trace 
=  magnification  at  period  T. 


(2-12) 


Ap,  is  easy  to  calculate  during  the  normal  procedure.  We 
shall  show  that  this  quantity  is  related  to  the  root  mean 
square  horizontal  amplitude  of  the  Rayleigh  waves: 

Assume  that  all  vertical  motion  is  due  to  Rayleigh 
waves  coming  from  a  single  direction.  Then  Love  waves  will 
not  contribute  to  A^  and 


A 


R 


0.798 


(2-13) 


will  be  an  estimate  of  the  mean  amplitude  of  the  envelope 
of  the  Rayleigh  wave.  This  could  yield  an  estimate  a  little 
too  large  since  a  small  crest  will  rarely  occur  among  the 
"well  developed"  crests  mentioned  in  the  procedure. 

The  time  interval  used  for  the  determination  of  the 
direction  of  approach  is  short,  and  hence  according  to  Longuet- 
Higgins  (Vesiac  Staff,  1962,  chapter  8)  the  microseisms  in 
general  and  the  Rayleigh  wave  in  particular  can  be  treated  as 
gaussian  noise,  and  hence(loc .  cit.)  the  quantity 
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height  of  crest  (from  mean  level) 

root  mean  square  amplitude 


(2.14) 


will  for  e  =  0,  i.e.  for  equal  number  of  zero  up  crossings 

and  crests,  follow  a  Rayleigh  distribution.  Hence  the  mean 
value  of  the  envelope  will  be 


and  now  it  follows  from  (2-13),  (2-l4)  and  (2-15)  that  AR 
will  be  an  estimate  of  the  root  mean  square  horizontal 
Rayleigh  amplitude. 


If  AR  is  small  compared  with  the  vertical  rms 


amplitude,  the  p  can  also  be  expected  to  be  small,  and  such 
an  observation  can  be  explained  by  assuming  that  a  substantial 
portion  of  the  vertical  amplitude  is  due  to  isotropic  micro¬ 
seisms;  an  alternate  explanation  is  that  Rayleigh  waves  are 
coming  in  from  two  directions  approximately  l80°  apart. 

2.5c  Reliability  of  small  p*s 

The  formula  (2-11)  gives  information  about  the 
probable  errors  as  determined  by  an  empirical  method  (Jensen, 
1958,  p.  8),  but  it  does  not  state  the  probability  distribu¬ 
tion  of  errors.  If  a  realistic  model  for  the  method  is 
used,  the  mathematical  problems  of  determining  this  prob¬ 
ability  distribution  are  very  involved. 
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A  simple  model  was  chosen,  however,  for  the 
purpose  of  determining  which  p’s  to  expect  from  a  random 
distribution  of  points.  For  the  model  21  points  are  consid¬ 
ered,  while  the  number  of  points  used  in  actual  practice 
often  exceeded  21  (see  Table  III). 

The  model  is  this:  21  Bernoulli  trials  are  made 
with  probability  for  "success".  If  the  number  of  "successes", 
s  ,  exceeds  10,  they  correspond  to  points  falling  in  the 
"right"  half  plane;  then 


I  lsn  ~  (21-sn)]  ‘  10 

I  21 


(2-16) 


will  correspond  to  the  p  defined  in  (2-10). 

If  the  number  of  "successes"  is  10  or  less  they 
correspond  to  points  in  the  "wrong"  half  plane  and  (2-l6) 
will  still  correspond  to  equation  (2-10)  . 

The  q’ s  defined  here  are  always  positive,  and  they 
are  generated  by  a  process  which  has  some  similarity  with  the 
positive  p’s  obtained  by  applying  Jensen’s  method  to  slopes 
measured  at  random  times,  rejecting  all  the  determinations 
with  p  ^  0. 

To  get  information  about  the  probability  distri¬ 
bution  of  q  we  will  apply  the  de  Moivre-Laplace  limit  theorem 
(Feller,  1950,  p.  137*)  which  for  the  present  purpose  yields 
a  good  approximation: 

pr  (a  ^  sn  £  6)  ~  $(xp+i/2)  "  $(xa-l/2^  (2"1^) 
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In  the  present  case 

h  =  (npq)-1/2  =  (21  •  •  ■|)'1//2  ~  0.436  (2-18) 

xt  =  (t-np)h  =  (t-10.5)  •  0.436  (2-19) 

Hence 

pr(q  ^  fy)  =  Pr(9^sn^12)  =  2[$(0. 872)-  |]  ~  0.61 

(2-20) 

or  q  >  1.4,  I.e.  q  ^  2.4  In  39 %  of  the  cases.  Similarly 
we  get 

q  >  3.3.,  I.e.  q  ^  4.3  In  8^  of  the  cases, 
q  >  5-2,  i.e.  q  ^  6.2  in  0.9^  of  the  cases. 

If  p  Is  equal  to  6  or  larger.  It  is  for  practical 
purposes  possible  to  conclude  that  the  noise  is  not  isotropic. 
If  two  determinations  with  p  ^  6  give  the  same  result  within 
the  limits  of  accuracy  it  is  possible  to  infer  with  complete 
confidence,  that  this  is  a  direction  of  approach,  provided 
that  the  basic  assumptions  of  the  method  are  satisfied. 

When  p  ^  3  no  conclusion  can  be  made  from  a  single  deter¬ 
mination,  but  the  distribution  of  such  directions  gives  a 
reliable  picture  of  the  most  frequent  directions  of  approach. 
As  indicated  by  the  probable  error  the  result  is  often  quite 
good  for  very  small  p’s,  which  shows  that  isotropic  noise 
is  not  a  common  phenomenon. 


, 
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Ill .  INSTRUMENTATION 

3 . 1  Short  Description  of  the  Instruments  Used 

The  data  used  in  the  preparation  of  the  present 
thesis  are  recorded  partly  by  the  University  of  Alberta 
observatory  instruments,  and  partly  by  the  portable  Texas 
Very  Low  Frequency  Refraction  System. 

The  University  of  Alberta  observatory  instruments 
were  installed  in  the  spring  of  1962  and  were  adjusted  to 
standard  frequency  response  and  calibrated  by  Mr.  F.  Lombardo, 
Dominion  Observatory,  using  a  Willmore  bridge.  The  instru¬ 
ments  consist  of  two  sets,  a  long  period  and  a  short  period 
set.  The  long  period  set  is  a  homogeneous,  three  component 
set  of  Press-Ewing  instruments,  manufactured  by  Lehner  and 
Griffith,  having  the  period  of  seismometers  adjusted  to  30 
seconds  and  the  period  of  the  galvanometers  to  90  seconds. 

The  short  period  set  consists  of  three  Willmore  seismometers, 
two  horizontal  Mark  I  and  one  vertical  Mark  II  used  with 
1/4  second  galvanometers.  All  the  components  are  recording 
on  Lehner  and  Griffith  triplet  recorders  which  usually 
eliminates  any  need  of  making  corrections  when  comparing 
two  components  of  a  set. 

Two  identical  sets  of  VLF  portable  seismic  record¬ 


ing  equipment  were  available.  One  of  these  sets  was  used 
by  the  author  for  studying  particle  motion  in  short-period 
microseisms.  The  equipment  was  installed  in  a  panel  truck. 


' 
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was  run  on  12  volt  batteries  and  It  consisted  of  a  100  volt 
power  supply,,  a  12  channel  amplifier  bank  with  adjustable 
attenuation  and  high  cut  filters,  a  RS-8U  camera  and  12 
vertical  seismometers  (2  cps) .  All  this  equipment  was  manu¬ 
factured  by  Texas  Instruments.  Two  Willmore  Mark  I  seismo¬ 
meters  (l  cps)  were  used  for  recording  horizontal  motion. 
Three  cables  with  4  take-outs  on  each  made  simultaneous 
recordingson  a  one  mile  spread  possible.  A  more  detailed 
discussion  of  the  field  equipment  is  given  in  the  manuals 
and  by  Weaver  (1962)  who  also  describes  some  of  the  problems 
that  arose  with  the  equipment  during  the  first  field  season. 

3.2  Determination  of  Phase  and  Amplitude  Response  of  the 
Field  Instruments 

During  the  field  season  1962  it  was  impossible  to 
obtain  a  homogeneous  set  of  instruments,  but  two  horizontal 
Willmores  with  periods  of  one  second  could  be  set  up  and 
compared  with  a  vertical  S-36  Texas  seismometer,  with  a 

•fcVie 

period  of  one-half  second.  Thus  a  knowledge  of,  phase  differ¬ 
ence  between  Willmore  and  Texas  seismometer  outputs  is 
necessary  in  order  to  apply  the  modified  Jensen  method.  In 
order  to  determine  this  phase  difference  for  simple  harmonic 
ground  motion,  let 


y 


yQ  +  a  sin  cot 


(3-D 
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be  the  coordinate  in  the  direction  of  recorded  motion  of  a 
point  of  the  frame  (ground) ,  let  x  be  the  coordinate  of  a 


point  of 

the  seismometer  mass  and  let 

(y  -  y0)  -  (x  -  xQ)  =  u 

(3-2) 

Then 

u  =  b  s  in  ( cot  -  6 ) 

(3-3) 

where  6  is  the  phase  angle  by  which  the  relative  motion  of 
the  ground  with  respect  to  the  seismometer  mass  lags  behind 
the  true  ground  motion. 

Let  the  force  acting  on  the  mass  M  be  Pu  +  Ru. 


Let 


P  2  R 

—  =  n  — 

M  M 


2A 


h 


Then  the  equation  of  motion  is 

u  +  2Au  +  n2u  =  y  (3-^0 

Now  from  (3-1)  (3-3)  and  (3-^): 

b(n2-0D2)  s in ( o)t - 6 )  +2Abca  cos(cot-6)  =  ao)2  sin  o)t  (3-5) 


For  the  three  particular  values  t  =  —  , 


O) 


and  t  =  0  the  equation  (3-5)  reduces  to 


■  a  2A  b 
Sln  5  =  —  •  - 


cos  6  = 


2  2 
n  -ao 


or 


b 

a 


2Aco 

2  2 
n  -a) 


t  = 


6 

CJD 


(3-6) 

(3-7) 


tan  6 


(3-8) 
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respectively.  Now  from  (3-6)  and  (3-7): 


(n2-cc>2)2  +  47^2^2 


(3-9) 


is  the  dynamic  magnification. 


In  the  Texas  S-36  geophone  a  15  K2  damping  resis¬ 


tance  is  built  in,  resulting  in  a  6K2  external  load.  This 
according  to  the  manual  corresponds  to  h  =  0.5  or  A  =  2tt. 
When  h  =  1  the  damping  is  critical,  and  h  =  0.5  corresponds 
to  a  damping  ratio  of  1:6.  With  these  particular  constants, 
calling  the  frequency  f,  equation  (3-8)  can  be  written: 


2f 


tan  6 


With  the  Willmore  Mark  I  seismometer  a  10002! 


damping  resistance  was  used.  The  average  damping  ratio  was 
measured  to  be  1:3  1/2,  corresponding  to  h  =  0.369. 


Using  these  constants  the  phase  difference  between 


the  output  of  the  two  seismometers  was  calculated  and  the 
results  are  shown  in  graphical  form  in  figure  4.  An  experi¬ 
mental  test  on  this  theoretical  curve  was  made  by  recording 
7  1/2  cps  microseisms  (at  University  campus,  night-time)  with 
a  vertical  Willmore  on  one  trace  and  a  Texas  S-36  seismo¬ 


meter  on  another  trace.  The  average  of  15  readings  was 
4.3  ±  1.2  milliseconds  corresponding  to  11°. 6  +  3°. 2,  with 
the  Willmore  output  lagging  behind  the  Texas  output  as 
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FREQUENCY  (CPS) 


S  :  Phase  laq  recorded  by  Willmorel  seismometer  (  I  cps  )  . 
w  3 


Sr  :  Phase  lag  recorded  by  Tex  as  (  S-36 )  seismometer  {  Zcps )  . 


FIG. 4. 


PERCENT  OF  ONE  CYCLE 
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predicted.  The  result  is  in  good  agreement  with  the  theo¬ 
retical  result,  which  is  10°. 4.  As  usual  in  geophysical 
literature  the  ±  sign  stands  for  the  estimated  standard 
deviation . 

The  next  logical  step  after  a  comparison  of  the 
phase  difference  between  two  outputs  is  to  compare  their 
amplitudes.  The  e.  m.  f.  induced  in  the  coil  of  a  seismo¬ 
meter  is  equal  to  a  constant  K  times  the  velocity  of  the 
coil  with  respect  to  the  magnet  so  with  the  harmonic  motion 
given  by  (3-3)  the  e.  m.  f.  is  Kbco  (mean  level  to  peak). 

When  coil  resistance  is  r  and  external  load  is  R  the  voltage 
U  across  the  outlets  of  the  seismometer  is 


RKbco 

R+r 


(3-10) 


The  dynamic  magnification  may  be  rewritten  from  (3-9) 


-  =  [(1-  AH2  +  4h2  £sf  l'1/2  (3-11) 
a  f  2  f d 

fn  is  the  natural  frequency  of  the  seismometer,  it  equals  7^. 
The  maximum  speed  of  ground  motion  is  called  v  (reached  twice 
in  a  cycle);  from  (3-1) 


v  =  aoo  =  2mfa 

Then,  combining  (3-10)  *  (3-H)  an<^  (3-12): 

R 


(3-12) 


U  =  K  * 


R+r 


[(i-in!)2  +  4h  2in!f1/2.v  0 


f‘ 


f‘ 
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This  formula  is  somewhat  similar  to  the  one  given  in  the  manual. 
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The  system  responds  linearly  to  the  voltage  U,  but 
the  amplification  depends  on  frequency,  high  cut  filter 
setting  and  the  attenuation  setting  used  on  the  amplifier. 

For  the  Texas  S-36  seismometer  the  constant  K  is 
given  in  the  manual  to  be  8.33  volts/inch/sec  or  3.28  volts^/cm/sec . 
The  total  external  load  for  the  Texas  seismometer  coil  when 
used  with  the  Very  Low  Frequency  Refraction  System  is  6KH, 
the  coil  resistance  is  40009. 

To  determine  the  value  of  K  for  the  Willmore  seismo¬ 
meter  a  comparison  of  Willmore  output  with  Texas  output  for 
20  cps  was  made.  First  one  sample  of  such  microseisms 
(Physics  Building,  basement)  was  recorded  with  Texas  seismo¬ 
meters  on  traces  1  and  2,  then  one  recording  was  made  with 
the  Texas  seismometer  (No.  1023)  replaced  by  Willmore  Mark  I 
(No.  84694),  set  up  for  vertical  operation.  By  comparisons 
with  the  common  seismometer  on  trace  2  it  was  possible  to 
make  the  comparison  independent  of  differences  in  amplification 
of  the  two  amplifiers.  The  result  was  that  K  for  the  Willmore 
I  seismometer  is  0.955  volts^/cm/sec .  For  high  frequencies, 
i.e.,  when  the  dynamic  magnification  is  1,  the  ratio  of  the 
Willmore  amplitude  to  the  Texas  amplitude  with  damping 
resistance  used  will  hence  be  0.313.  For  lower  frequencies 
the  difference  between  the  frequency  responses  of  the  two 
instruments  must  be  taken  into  account.  The  change  of  trace 
amplitude  caused  by  this  effect  is  shown  on  figure  5;  the 
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Normalized  velocity  sensitivity  of  Texas  S-36  seismometer 


Normalized  velocity  sensitivity  of  Willmore  I  seismometer 


Normalized  ratio 


Texas 

Will  more 


of  velocity  sensitivities 


Normalized  d:  Equals  one  at  high  frequencies 


FIG. 5. 
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response  is  normalized  to  be  one  for  high  frequencies. 

The  input  voltage  to  amplifiers  can  be  derived  from 
the  trace  amplitude  by  the  formula: 


log  U  = 


Si 

log  A  +  — £  +  log  B  +  log  C 


(3-14) 


Where  U  -  input  to  amplifiers  in  volts  (measured  from  peak 
to  average  voltage) 

A  -  trace  amplitude  in  millimeters  (peak  to  peak) 

B  -  correction  factor  due  to  insertion  of  filters 
C  -  a  constant  of  the  amplifier  system,  implicitly 

_7 

given  in  manual  to  be  2.95  *  10 
a,p  -  attenuation  setting,  variable  from  0  to  66  db  of 
energy  attenuation. 

Table  1  shows  log  B  as  a  function  of  filter  settings  and  of 
the  frequency;  the  values  are  derived  from  the  filter  response 
curve  in  the  manual.  Combining  (3-11) >  (3-13)  and  (3-14) 
yields : 


log  v  =  log  A  +  —  +  log  B  +  log[d.|±d  .  _1_]  (3-15) 


This  formula  shows  that  the  dynamic  magnification  b/a  is  a 
measure  of  the  velocity  sensitivity.  Since  b/a  equals  one 
at  high  frequencies,  it  can  be  called  the  normalized  velocity 
sensitivity.  The  values  of  this  quantity  are  shown  in  graph¬ 
ical  form  in  figure  5.  Combining  (3-12)  and  (3-15): 


=  log 


a  rn 

A  +  — 
20 


+  log  B  +  log[ 


C(R  +  r)  j 

KR  •  2tt  •  b/a 


log  f 

(3-16) 


log  a 
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The  formula  (3-15)  for  the  ground  velocity  Is  some¬ 
what  simpler  than  (3-16)  Inasmuch  as  the  frequency  dependence 
only  shows  In  the  factor  b/a  as  a  low-cut  filter  effect  and 
in  the  factor  B  as  a  high-cut  filter  effect;  in  between  the 
response  is  flat.  Thus  it  is  easier  to  compare  amplitudes 
connected  with  different  frequencies  when  ground  velocities 
are  used,  and  also  the  smallest  measurable  seismic  noise  will 
have  a  flat  frequency  dependence  in  the  range  considered. 

The  use  of  ground  velocities  to  measure  vibrations  can  also 
be  considered  as  a  compromise  between  the  use  of  ground 
amplitude  ( oc  ground  velocity  times  f”1)  and  ground  accelera¬ 
tions  (oc  ground  velocity  times  f)  . 

Measuring  a  in  Angstroms  and  v  in  Angstroms/sec, 
the  numerical  value  of  K  will  be  modified  so  that  the  factor 


attains  the  values:  15.0  for  the  Texas  seismometer 
and  4.70  for  the  Willmore  seismometer. 

The  resulting  values  of  the  quantity  log  -L  • 


1 

b/a 


are  shown  in  Table  H  for  the  Texas  S-36  seismometer. 

To  get  the  ground  velocity  in  Angstroms/sec, 
formula  (3-15)  is  convenient  to  use  with  Tables  1  and  2. 
Ground  amplitude  if  desired  may  then  be  obtained  by  (3-12)  or 


(3-16) . 

The  magnifications  derived  from  the  above  formulas 
are  only  half  of  those  obtained  by  means  of  Dr.  Garlandfe 

(Weaver  1962) .  The  reason  is  that  in  this  work  the  use 


curve 
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of  peak  to  peak  amplitudes  was  preferred  for  the  trace 
amplitudes,  hut  not  for  the  velocities  of  the  ground  motion. 


46. 


TABLE  I 

Log  Bj  the  Reduction  for  the  Attenuation  due  to  Filters. 


T 


‘Frequency 

CPS 

12 

—  r lulu 

16 

24 

32 

48 

2  | 

0.26 

0.22 

0.17 

0.10  ' 

0.05 

0.00 

3 

0.28 

0.22 

0.17 

0.10 

0.05 

0.00 

4 

1  0.31 

0.23 

0.17 

0.10 

0.05 

0.00 

5 

:  0.36 

0.25 

0.17 

0.10 

0.05 

0.00 

6 

0 .4l 

0.28 

0.17 

0.10 

0.05 

0.00 

7 

j  0.47 

| 

0.31 

0.19 

0.11 

0.05 

0.00 

8 

0.53 

0.35 

o.2l 

0.12 

.  .  .  ... 

0.06 

0.00 

9 

0.57 

0.38 

0.23 

0.13 

0.06 

0.00 

10 

o .  66 

0.42 

0.25 

0 . 14 

0.06 

0.01 

12 

0.78 

0.49 

0.30 

0.16 

0.07 

0.01 

14 

0.92 

0.58 

0.36 

0.19 

0.08 

0.02 

16 

1.66 

0 . 66 

0 .41 

0.22 

0.10 

0 .04 

18 

1.19 

0.74 

0 .46 

0.26 

0.12 

0.05 

20 

1.31 

0.82 

0.52 

0.29 

0 . 14 

0.06 

22 

0.90 

0.58 

0.33 

0.16 

0.08 

24 

0.98 

0.63 

0.37 

0.19 

0.09 

26 

1.05 

0.69 

0 . 4l 

0.21 

0.11 

28 

1.12 

0.74 

0.45 

0.24 

0.13 

30 

1.19 

0.80 

0.50 

0.27 

0 . 14 

35 

1.37 

0.95 

0.59 

0.33 

0.19 

4o 

1.55 

1.08 

0.70 

0.40 

0.24 

45 

1.21 

0.8l 

0.48 

0.28 

50 

1 

J 

0.90 

0.55 

0.34 

t 
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TABLE  II 


Values  of  the  Quantity  log  —  for  the  Texas  S-36 

KR  •  b/a 

Seismometer,  to  be  used  with  Formula  (3-15). 


Frequency 

1 

1.5 

2 

2.5 

— 

3 

4 

5 

— 

6 

- - 1 

7 

log  cjR±h_ 

KR  •  b/a 

1.7^ 

1.37 

1.18 

C\J 

1 — 1 

• 

1 — 1 

1.12 

-  

1.14 

1.15 

1.16 

1.16 

Frequency 

8 

9 

10 

11 

12 

13 

14 

>14 

loghR±h_ 

KR  •  b/a 

1.16 

1.17 

1.17 

1.17 

1.17 

1.17 

1.18 

1.18 
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IV.  THREE  TO  TEN  SECOND  MICROSEISMS  AT 

STATION  EDMONTON  JULY  AND  AUGUST  1963 

4 . 1  Introduction 

Most  microseisms  in  the  3  to  10  second  period 
range  are  caused  by  disturbed  weather  over  large  bodies  of 
water.  They  show  a  clear  seasonal  variation  resulting  in 
larger  mean  periods  and  mean  amplitudes  in  winter  as 
compared  to  the  summer.  The  distribution  of  amplitudes 
in  a  given  month  is  markedly  skew,  making  it  possible  to 
distinguish  between  microseismic  storms  and  the  usual  back¬ 
ground  microseisms.  The  term  "microseismic  storm"  is 
used  to  describe  a  period  where  the  amplitudes  are  larger 
than  normal . 

Discussion  on  various  general  aspects  of  the  3  to 
10  second  microseisms  was  given  in  the  introduction,,  and 
further  discussions  will  follow  in  section  4.5a.  The 
present  chapter  will  chiefly  be  devoted  to  the  inter¬ 
pretation  of  microseisms  recorded  at  station  Edmonton  (the 
University  of  Alberta  seismic  observatory) . 

4.2  Amplitudes  and  Periods 

For  the  two  months,  July  and  August  1963,  the 
amplitudes  of  3  -  6  second  microseisms  at  station  Edmonton 


t 
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on  the  LPZ  (i.e.  long  period  vertical  seismograph)  were 
measured  at  6  hour  Intervals .  Each  amplitude  was  the  aver¬ 
age  of  maximum  amplitudes  in  ten  minute  intervals.  The 
results  are  given  in  figure  6.  Most  of  the  periods  were 
between  4  and  5  seconds  and  therefore  it  was  found  reasonably 
accurate  to  use  a  common  factor  to  convert  peak  to  peak  trace 
amplitude  to  center  to  peak  ground  amplitude.  The  factor 
used  corresponded  to  the  period  of  4  1/2  seconds. 

It  will  be  seen  from  figure  6  that  it  is  easy  to 
locate  microseismic  storms;  the  most  conspicious  ones 
occurred  at  July  6  and  at  August  13,  15  and  28. 

Only  on  two  occasions  did  periods  longer  than  6 
seconds  occur  with  an  amplitude  exceeding  that  of  the 
simultaneous  3-6  second  microseisms.  One  occasion  was 
August  6  with  a  few  beats  of  8  second  microseisms  having  an 
amplitude  of  0.l4  microns.  The  other  occasion  was  a  micro- 
seismic  storm  of  8  second  microseisms .  This  storm  was 
not.icable  from  August  23,  9h  GMT  to  August  26,  lh  GMT.  The 
amplitude  was  gradually  increasing  through  the  value  0.21 
microns  at  August  24,  6h  GMT  to  0.35  microns  at  August  24, 
l8h  GMT.  These  8  second  microseisms  are  not  shown  in 

figure  6. 

The  period  of  the  dominant  microseisms,  i.e.  of 
the  microseisms  that  produced  well  developed  maxima  on  Z, 
were  measured  whenever  a  direction  of  approach  was  determined 


•i 
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FIG.  6.  Amplitudes  of  3-6  second  microseisms  on  the  longperiod  vertical  (  LPZ)  recorded  at  stati 
Edmonton, the  University  of  Alberta  seismic  observatory ,  during  July -August  1963. 
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(see  section  4.4).  The  results  of  these  period  measurements 
are  shown  in  Table  III. 

4 . 3  The  Weather  Situation  in  July  and  August  1963 

A  complete  description  of  what  is  known  about  the 
weather  situation  along  the  North  American  coasts  during 
July  and  August  1963  would  take  several  books .  It  is  avail¬ 
able  on  punched  cards.  Some  general  features  of  the  weather 
situation  associated  with  periods  of  anomalous  activity  of 
microseisms,  however,  could  be  extracted  by  looking  through 
the  weather  maps.  These  maps  were  kindly  made  available  to 
me  by  E.  Reinelt,  Meteorological  Service  of  Canada,  Edmonton 
Municipal  Airport.  The  information  given  in  this  section 
will  be  used  in  the  discussions  in  section  4.5.  The  winds 
quoted  are  winds  reported  from  weather  stations,  generally 
the  highest  wind  reported  from  the  area  in  question. 

July  5,  6h  GMT .  A  998  mb  barometric  low  is 
situated  near  Pt .  Barrow  moving  ENE. 

July  6,  6h  GMT.  The  low  has  deepened  to  998  mb, 
and  is  now  at  75°N,  121°W  (McClure  Strait) .  For  the  season 
it  is  quite  deep,  and  it  is  producing  30  knot  winds.  The 
other  coastal  areas  are  fairly  quiet  during  this  period. 
Winds:  0-5  knots  at  the  Californian,  Oregon,  Washington 
and  B.C.  coasts;  15  knots  at  the  Alaskan  SW  coast. 


-  - '• 
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July  6,  18^  GMT .  The  low  has  moved  to  77°N, 

120°¥ . 

h. 

July  7^  6  GMT.  The  low  has  deepened  to  992  mb 
and  Is  now  situated  at  79°N,  110°¥. 

July  8,  0h  GMT .  The  low  has  weakened  to  1000  mb 
and  moved  to  79°^  120°¥. 

Tuly_13^  6h  GMT .  A  998  mb  low  in  the  Hudson  Bay 
is  producing  20  knot  winds. 

July  21,  l8h  GMT .  A  1004  mb  low  is  situated  at 
48°N,  132°¥,  4  degrees  off  the  coast  of  B.  C. 

July  22,  0h  GMT .  The  low  is  just  one  degree  off 
the  coast  and  has  deepened  to  1000  mb  and  is  producing  30 
knot  winds.  Maximum  wind  observed  in  Hudson  Bay,  The 
Canadian  Arctic  and  Labrador  is  15  knots. 

h 

July  22,  12  GMT.  The  low  has  moved  north  along 
the  coast  to  53°N  and  has  deepened  to  998  mb.  The  maximum 
wind  along  northern  and  eastern  Canadian  coasts  is  still 
15  knots. 

July  23,  _0^  GMT.  The  low  has  filled  up,  but  30 

knot  winds  still  occur. 

"h 

July  23,  6  GMT .  20  knot  winds  occur  at  the 


Pacific  coast  and  in  the  Hudson  Bay. 
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August  13,  6  GMT .  A  low  In  the  northern  part  of 
Labrador  Is  associated  with  25  knot  winds.  Along  the  coast 
of  B.  C.  20-40  knot  winds  occur. 

August  13,  l8  GMT .  A  frontal  system  now  developed 
in  the  Hudson  Bay  area  and  a  low  at  the  B.  C.  coast  are  both 
producing  20  knot  winds. 

August  l4,  12^  GMT.  A  1012  mb  low  at  68ON,  110°W 
(Coronation  Gulf)  producing  15  knot  winds. 

August  14,  l8^  GMT .  The  low  has  moved  to  68°N, 
106°W  and  winds  have  increased  to  25  knots.  The  winds  along 
the  B.  C.  coast  and  over  Hudson  Bay  are  calm  or  fairly  light. 

August  15,  6h  GMT .  The  low  is  1008  mb  at  72°N, 

95°W  with  25  knot  winds.  Another  low  is  situated  at  the 
northern  part  of  Labrador. 

August  15,  l8h  GMT .  The  Labrador  low  is  now 
producing  30  knot  winds . 

August  22,  0h  GMT .  Two  deep  lows  (both  996  mb, 
fairly  deep  for  a  summer  low)  are  associated  with  25-30 
knot  winds.  One  is  at  48°N,  l66°W  (the  Pacific,  5  degrees 
from  the  Aleutian  Islands)  slowly  moving  east,  while  the 
other  is  at  58°N,  l65°W  (the  Bering  Sea)  rapidly  moving 
towards  the  Bering  Strait. 
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August  22,  12*1  GMT.  The  lows  have  deepened  to 
992  mb  at  47°N,  l62°W  and  984  mb  at  62°N,  l67°W. 

h 

August  23j  0  GMT.  The  cyclones  are  still  very 
active;  the  Pacific  one  Is  1000  mb  at  50°N,  l6o°W  with  strong 
winds  (35  knots)  In  an  annular  pattern  3-5  degrees  from  the 
center  of  the  cyclone.  The  other  low  Is  now  984  mb  at  69°N, 
l67°¥  (north  of  the  Bering  Strait)  and  Is  producing  30  knot 
winds . 

August  23j  l8^  GMT .  The  Pacific  low  is  now  1004  mb 
and  at  53°N,  157 °¥  with  fronts  extending  10  degrees  S  and  SE 
of  that  point  and  with  25  knot  winds  at  these  fronts. 

August  24,  12h  GMT .  The  low  has  been  filled  up, 
but  20  knot  winds  still  occur  in  the  Pacific. 

August  28,  0h  GMT .  A  low  in  Baffin  Bay  produces 
25  knot  winds  in  Greenland.  The  Canadian  Arctic  is  also 
disturbed  with  25  knot  winds  at  Coronation  Gulf. 

August  28,  6h  GMT .  The  Baffin  Bay  low  produces 
25  knot  winds  in  Greenland,  a  front  system  in  the  Canadian 
Arctic  produces  25  knot  winds  at  120°W . 

August  28,  12h  GMT.  The  low  in  Baffins  Bay  pro¬ 
duces  25  knot  winds  in  the  straits  north  of  Hudson  Bay. 

The  Arctic  front  system  produces  20  knot  winds  at  Coronation 


Gulf. 
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4 . 4  The  Directions  of  Approach 

The  directions  of  approach  were  determined  in  2 6 
cases  selected  from  the  months  July  and  August  1963  at 
station  Edmonton.  The  method  used  was  Jensen’s  method., 
described  in  section  2.5.  Altogether  1258  slopes  were  read. 

The  results  are  shown  in  Table  III,  and  some  of 
the  NE  diagrams  used  during  the  determinations  are  shown  in 
Figure  7.  T  is  the  period  as  measured  during  the  time 
interval  used  for  the  determination  of  the  direction  of 
approach;  a  is  the  direction  of  approach,  i.e.  the  directions 
that  the  waves  come  from,  measured  from  north,  positive 
towards  east.  The  reliability  measure  p  and  the  "Rayleigh 
amplitude"  AR  are  defined  in  section  2.5.  The  significance 
of  the  three  different  categories.  A,  B,  and  C  is  given 
below;  n  is  the  number  of  points. 

At  first  11  determinations  were  made  in  order  to 
find  the  directions  of  approach  of  the  microseismic  storms. 
Two  or  three  determinations  were  made  on  each  storm,  and  the 
results  agreed,  so  that  the  directions  determined  for  a 
given  storm  always  belonged  to  the  same  or  to  two  adjacent 
octants.  The  12  determinations  of  the  approach  of  storm 
microseisms  are  labeled  "category  A". 

The  weather  maps  were  studied,  and  the  correlation 
with  wind  storms  in  the  Arctic,  in  Hudson  Bay  and  at  Labrador 
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TABLE  III 


The  Directions  of  Approach  of  Microseisms  at  Edmonton 
Determined  by  Jensen’s  Method. 


DATE 

Time  Inter¬ 
val  used 
h  m 

T 

(sec) 

p 

a 

cate¬ 

gory 

n 

ar 

(milli¬ 

microns) 

Oc¬ 

tant 

July  2 

1 

10  1 

01-05 

4.3 

4 

+159° 

C 

20 

18 

S 

6 

5  ; 

02-08 

5.1 

10 

+  12° 

A 

19 

288 

N 

8  ! 

53-56 

5.1 

9 

+  12° 

A 

21 

244 

N 

18 

01-07 

4.8 

9 

+  31° 

A 

20 

172 

NE 

7 

6  1 

01-08 

5.2 

10 

+  54° 

B 

23 

91 

NE 

10 

10 

01-06 

3.9 

4 

0° 

C 

26 

14 

N 

12 

10 

01-05 

4.1 

3 

+110° 

C 

25 

20 

E 

13 

6 

01-06 

4.0 

5 

+108° 

B 

35 

26 

E 

6 

06-10 

3.9 

1 

+148° 

B 

26 

18 

SE 

19 

10 

01-06 

4.1 

1 

-166° 

C 

29 

7 

S 

22 

9 

01-04 

4.0 

2 

+  44° 

A 

24 

30 

NE 

11 

56-63 

4.4 

0 

+  90° 

A 

24 

37 

E 

23 

10 

01-05 

4.0 

5 

+  63° 

C 

30 

33 

NE 

24 

10 

01-06 

3.5 

4 

+  37° 

C 

27 

22 

NE 

29 

10 

01-07 

4.1 

3 

+  56° 

C 

24 

14 

NE  1 

1 

Aug.  3 

10 

01-06 

4.6 

1 

+  97° 

c 

28 

7 

E 

4 

10 

01-07 

4.2 

5 

+  78° 

c 

26 

14 

E 

10 

10 

01-07 

3.8 

3 

+  42° 

c 

33 

11 

NE 

13 

6 

01-03 

4.4 

7 

+  68° 

A 

19 

95 

E 

9 

01-03 

4.5 

7 

+105° 

A 

1 

18 

128 

E 

12 

01-04 

4.4 

10 

+  50° 

A 

19 

108 

NE 

14 

17 

57-63 

4.5 

3 

+  72° 

1 

B 

18 

82 

E 

24 

5 

55-61 

\ 

1 

6 

27-28 

V  8.2 

9 

-  71° 

A 

24 

209 

¥ 

16 

53-59 

7.9 

9 

-  79° 

A 

1 

27 

164 

¥ 

28 

'  3 

01-04 

5.2 
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FIG.  7.  N E  -  diagrams, used  to  get  some  of  the  directions  listed  in 

table  m. 
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was  found  to  be  so  good  that  an  attempt  to  predict  directions 

of  approach  from  some  weather  situations  was  undertaken. 

These  directions  are  labeled  "category  B" . 

Finally  ten  dates  were  selected  among  the  remaining 

dates ,  i.e.  among  the  quiet  periods  between  the  microseismic 

storms.  A  table  of  random  numbers  (Hald,  1951)  was  used  to 

select  these  dates.  The  record  amplitude  at  these  dates  was 

only  a  small  fraction  of  a  millimeter.  The  directions  of 

h 

approach  at  these  dates  was  determined  at  the  time  10  GMT, 
where  the  long  period  noise,  presumably  of  thermal  origin, 
was  a  low  level.  These  ten  directions  are  labeled  "category 
C". 


4 . 4a  Distribution  into  Octants 

Table  IV  shows  the  distribution  into  octants  for 
observations  having  different  period  ranges.  Category  C 
observations  are  listed  separately. 

The  directions  of  approach  of  the  3-5  second 
microseisms  scatter  chiefly  in  the  octants  E  and  NE,  and 
this  tendency  is  the  same  for  the  microseismic  storms  and 
for  the  low  activity  microseisms. 

The  5-6  second  microseisms  came  from  N  and  NE 

and  all  the  data  originate  from  two  storms. 
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TABLE  IV 

Joint  Distribution  Period  T  in  seconds  versus 
the  Octant  of  the  Direction  of  Approach. 


CATEGORY 

C 

CATEGORY  A  and  B 

TOTAL 

3<T<5 

3<T<5 

5<t<6 

6<T<7 

7<T<9 

N 

1 

0 

3 

0 

0 

4 

NE 

4 

3 

2 

0 

0 

9 

E 

3 

5 

0 

0 

0 

8 

SE 

0 

1 

0 

0 

0 

1 

S 

2 

0 

0 

0 

0 

2 

sw 

0 

0 

0 

0 

0 

0 

w 

0 

0 

0 

0 

2 

2 

NW 

0 

0 

0 

0 

0 

0 

TOTAL 

10 

9 

5 

0 

1 _ 

2 
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No  6  -  7  second  .microseisms  were  observed  during 
the  two  months  under  study. 

The  7-9  second  microseisms  were  observed  with 
sufficient  amplitude  only  during  one  8  second  microseismic 
storm.  The  reliability  of  these  directions  is  very  good 
(p  =  9  in  both  determinations)  and  they  are  the  only  obser¬ 
vations  in  the  western  octant. 

4 . 5  Interpretation  of  the  Directions  of  Approach 

During  the  interpretation,  some  quantitative  work 
was  done  based  on  Longuet-Higgins'  theory;  however  only  in 
cases  where  this  theory  is  mentioned  does  the  conclusion 
depend  on  this  theory;  conclusions  concerning  distributions 
of  directions  of  approach  and  selective  barriers  are 
independent  of  assumptions  concerning  the  mechanism  of  the 
generation  of  microseisms. 

4 . 5a  Microseisms  Expected  from  Given  Wind  Velocities 
According  to  Longuet-Higgins*  Theory 

If  Longuet-Higgins1  theory  is  used  as  a  working 
hypothesis,  generation  of  3-8  second  microseisms  will  only 
be  expected  from  large  bodies  of  open  water,  and  this 
consequence  agrees  with  the  observations. 

Longuet-Higgins'  theory  was  further  developed  by 
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Hasselmann  (1963)^  who  took  the  areal  extent  of  the  source 
Into  account.  He  also  pointed,  out  that  microseisms  would 
he  produced,  not  by  interaction  between  gravity  waves  having 
exactly  equal  and  opposite  wave  numbers,  but  rather  by 
interactions  between  waves  with  slightly  different  wave 
numbers,  that  produce  compressional  waves  of  finite  hori¬ 
zontal  phase  velocity  in  resonance  with  a  trapped  mode  of 
the  waveguide . 

The  theory  predicts  proportionality  between  the 
amplitude  of  microseisms  and  the  square  of  the  waveheights 
of  the  interfering  waves.  This  prediction  was  confirmed 
by  Hjortenberg  (1963).  The  square  can  be  taken  to  be 
proportional  to  the  energy  of  the  waves,  and  this  energy, 
according  to  Defant  (1961,  p.  55)  is,  for  a  fully  arisen 
sea,  proportional  to  the  fifth  power  of  the  wind  velocity. 

The  periods  of  wind  waves  to  expect  from  a  given 
wind  velocity  of  long  duration  and  large  horizontal  extent 
was  given  implicitly  by  Defant  (1961,  p.  5*0  • 

T  =  0.40  v  (4-1) 

max 

where  T  is  the  water-wave  period  associated  with  the 
max 

peak  of  the  power  spectrum  and  v  is  the  wind  speed  in 
knots  . 

Hence  20  and  30  knot  winds  would  produce  water 
spectra  with  peaks  at  8  and  12  seconds,  and  from 


wave 
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Longuet-Higgins ’  theory  4  and  6  second  period  .microseisms 
would  he  expected.  (Longuet-Higgins,,  1950). 

This  prediction  agrees  with  the  microseismic 
storms  investigated  in  Edmonton.  The  periods  associated 
with  category  A  observations  of  non-Pacific  origin  were 
4.0  -  5.8  seconds  and  the  winds  observed  in  the  suspected 
source  areas  were  20  -  30  knots. 

4 . 5^  Comparison  of  Microseisms  with  High  Winds  in  the 
Pacific 

Three  cases  of  disturbed  weather  in  the  Pacific 
are  described  in  section  4.3.  The  first  case  was  between 
July  22,  0h  GMT  and  July  23.,  0h  GMT.  Winds  were  30  knots, 
while  winds  elsewhere  were  less  than  15  knots;  it  should 
be  mentioned,  however,  that  the  weather  maps  available  did 
not  permit  investigation  of  the  weather  along  the  eastern 
coast  of  U.  S.  A.  The  directions  of  approach  determined 
are  in  the  octants  E  and  NE,  but  the  pls  are  very  small 
and  figure  7  shows  that  the  points  are  scattered  in  all 
directions.  The  AR  is  only  one  third  of  the  AR  determined 
on  August  13,  at  which  time  the  vertical  amplitude  was 
similar.  This  could  indicate  that  microseisms  of  Pacific 
origin  account  for  roughly  one-third  of  the  amplitude. 

From  the  weather  situation,  however,  the  Pacific  microseisms 
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could  be  expected  to  exceed  other  microseisms.  From 

Longuet-Higgins 1  theory  they  could  be  expected  to  have  at 
r 

least  2  =  32  times  the  amplitude  from  other  sources,, 

but  their  amplitude  was  smaller,  and  this  Is  Interpreted 
as  a  very  heavy  damping  of  Pacific  microseisms. 

Yi 

The  second  case  is  August  13,  6  GMT.  This  was 
a  storm  affecting  local  areas  of  the  Pacific  coast  and  was 
not  associated  with  any  extended  low  pressure  system.  It 
appears  from  figure  7  that  no  perturbation  of  the  .micro- 
seisms  were  produced. 

The  third  case  at  August  22  to  August  23  was 
an  extended  low  pressure  system  in  the  Pacific.  It  will 
be  discussed  in  the  next  subsection. 

4.5c  Eight  Second  Microseisms 

Only  one  microseismic  storm  of  8  second  period 

occurred  during  ih  e  two  months  investigated.  The  amplitudes 

h. 

are  described  in  section  4.2.  From  August  24,  6  GMT  to 
August  24,  l8h  GMT,  these  microseisms  dominated  over  shorter 
period  microseisms,  and  the  directions  of  approach  could 
be  determined.  The  amplitude  was  steadily  increasing  until 
August  24,  l8h  GMT.  For  the  next  24  hours  no  record  was 
available  and  after  that  the  amplitudes  were  again  very 
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Both  directions  of  approach  were  characterized  by 
p  =  9  and  they  were  pointing  towards  the  B.  C.  -  Alaska 
boundary  at  the  Pacific  coast. 

The  weather  situation  in  the  Pacific  had  no 
unusual  features  during  the  time  of  the  8-second  microseismic 
storm,  so  the  weather  during  this  period  could  not  explain 
it. 

Now  let  us  apply  Longuet-Higgins *  theory  as  a 
working  hypothesis: 

Since  no  low  pressure  with  a  simultaneous  genera¬ 
tion  of  water  waves  with  different  directions  is  available, 
it  is  reasonable  to  attribute  the  microseisms  to  interfer¬ 
ence  between  in-coming  and  reflected  swells  with  a  period 
of  l6  seconds. 

For  most  of  the  path  the  water  is  deeper  than 
half  a  wavelength  and  the  amplitude  is  small  compared  to 
the  wavelength,  hence  from  the  work  of  Barber  and  Tucker 

(1962): 

V  =  (4-2) 

group  47 r 

Thus  the  group  velocity  of  a  l6  second  gravity 

wave  is  0.405  degrees  per  hour. 

According  to  Defant  (1961,  p.  53)  the  power 
density  of  the  water  wave  energy  at  a  period  of  16  seconds 
can  be  expected  to  be  8  1/2  times  the  power  density  when  the 

wind  velocity  is  40  knots  as  compared  to  the  power  density  when 
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the  winds  have  a  speed  of  30  knots.  During  August  22, 

6  -  18  GMT,  the  strongest  winds  reported  were  30  knots, 

but  at  August  23,  O^1  GMT  the  maximum  of  all  the  winds 
reported  from  the  low  occurred.  This  wind  speed  was  35 
knots,  and  an  annular  pattern  of  equidistant  isobars  made 
it  reasonable  to  assume  that  35  knot  winds  occurred  through¬ 
out  this  annular  pattern,  3-5  degrees  from  the  center. 

The  distance  from  the  northernmost  point  of  the 
B.  C.  Coast  to  the  closest  point  of  this  annular  pattern  is 
12  degrees.  This  corresponds  to  30  hours  of  travel  time 
for  the  16  second  water  waves,  so  these  waves  would  arrive 
at  August  24,  6h  GMT. 

The  distance  at  which  the  largest  amount  of  wind 
energy  in  the  annular  zone  has  a  direction  towards  the  con¬ 
tinent  will  correspond  to  the  distance  to  the  center  of  the 
low.  This  distance  is  17  degrees,  and  hence  these  large 
amplitude  waves  are  expected  to  arrive  42  hours  later,  i.e. 
at  August  24,  18^  GMT. 

h 

Now  the  interval  between  August  24,  6  GMT  and 
August  24,  l8h  GMT  is  exactly  the  interval  where  the  microseism 
amplitude  is  steadily  increasing  to  the  high  level  at  Edmonton. 

Furthermore,  35  knot  winds  would  according  to  (4-1) 
produce  water  waves  with  a  power  density  peak  at  14  seconds, 
and  according  to  Defant  (1961,  p.  53)  the  range  between  the 
power  points  is  approximately  one  octave,  so  the  16 
second  water  waves  are  well  within  the  period  range  expected. 
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The  consequences  from  the  use  of  Longuet-Higgins ’ 
theory  as  a  working  hypothesis  are  thus  in  good  agreement 
with  these  observations. 

A  swell  having  a  period  of  8  seconds  from  the  same 
source  cannot  account  for  the  observations,  because  it 
travels  at  half  the  speed  and  would  arrive  at  least  30  hours 
late.  It  is  also  impossible  that  the  same  low  could  have 
generated  the  8  second  water  waves  at  an  earlier  time, 
because  the  low  was  moving  toward  the  continent  with  a 
speed  approximately  equal  to  the  group  velocity  of  8  second 
gravity  waves. 

Since  wavegauge  recordings  were  not  available,  it 
is  impossible  to  prove  that  the  period  of  the  swell  was  l6 
seconds  and  it  can  only  be  stated  that  such  a  swell  is  to 
be  expected  from  the  weather  situation. 

One  additional  feature  of  the  8  second  storm 
may  be  worthwhile  to  point  out:  The  azimuth  at  station 
Edmonton  of  the  position  of  the  Pacific  low  at  August  23, 

0h  GMT  was  a  =  -  73°.  This  value  agrees  with  the  directions 
of  approach  of  8  second  microseisms  at  station  Edmonton  30 
to  4l  hours  later. 

This  agreement  can  be  explained  by  the  action  of 
swell  at  the  coast  using  the  modification  of  Longuet-Higgins 3 
theory,  introduced  by  Hasselmann  (1963)  3.nd  described  in 
section  4.5 s.,  and  one  additional  assumption. 
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The  argument  Is  this:  The  wavefronts  are  perpen¬ 
dicular  to  the  direction  of  wave  propagation,  and  waves 
reflected  from  the  coast  may  he  assumed  to  have  the  most 
constructive  interference  with  the  incoming  waves,  when 
the  direction  of  propagation  o f  the  reflected  waves  differs 
exactly  l80°  from  the  direction  of  the  incoming  waves.  This 
assumption  seems  reasonable  because  then  the  particle  motion 
is  in  phase  along  the  crest;  however,  it  is  nothing  but  an 
assumption . 

When  such  a  l8o°  difference  occurs,  all  the 
vectorial  wavenumbers  will  lie  on  a  straight  line  towards 
the  course,  and  so  will  the  differences  between  wavenumbers 
of  incoming  and  reflected  waves.  The  distribution  of  the 
absolute  values  of  wavenumbers  of  incoming  and  reflected 
waves  can  be  expected  to  be  identical,  and  thus  the  differ¬ 
ences  between  the  vectorial  wavenumbers  will  have  a 
symmetrical  distribution  on  the  line  towards  the  source  of 
the  swell.  A  certain  range  of  these  vectorial  wavenumbers 
will  satisfy  the  condition  to  set  up  seismic  waves,  and 
hence  the  seismic  energy  will  be  sent  out  in  two  directions; 
towards  the  source  of  the  swell  and  away  from  this  source. 

The  mechanism  described  here  can  explain  the 
direction  observed  in  Edmonton,  and  also  the  high  accuracy 
of  the  determination;  such  an  accuracy  would  result  from 
a  narrow  sector  of  approach. 
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A  fair  proportion  of  the  rocky  shores  of  the 
islands  in  the  suspected  source  area  for  the  microseisms., 
have  a  direction  perpendicular  to  the  direction  to  the 
source  of  the  postulated  swell. 

For  this  eight  second  storm,  Longuet-Higgins 1 
theory  has  explained  the  following:  The  period,  the 
variation  of  the  amplitudes  with  time,  the  direction  of 
approach  and  the  variability  of  the  same. 

The  evidence  given  here  favors  Longuet-Higgin ’ s 
theory  of  microseism  generation  for  this  particular  storm. 

4 . 5d  Comparison  Between  the  Directions  of  Approach  and 
the  Non- Pacific  Storms. 

The  Arctic  low  of  July  6  was  moving  along  the 
western  edge  of  the  Arctic  archipelago.  The  direction  to 
the  center  of  the  low  varied  between  a  =  -5°  and  cl  =  -1°, 
but  the  directions  of  approach  of  microseisms  in  Edmonton 
varied  between  +12°  and  +  31°. 

This  difference  may  be  explained  by  the  distribu¬ 
tion  of  the  polar  ice.  In  most  of  the  Arctic  Ocean 

the  water  is  covered  with  ice  all  summer,  and  the  winds 
would  not  produce  microseisms.  Between  the  islands  of  the 
archipelago  the  water  is  open  in  summer,  and  the  bodies  of 
water  are  large  enough  for  the  generation  of  microseisms. 
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At  July  7,  6  GMT,  the  low  was  deepening,  but  the  microseismic 
storm  had  a  small  and  decreasing  amplitude.  This  may  be 
due  to  the  position  ol  the  low  at  79°N;  at  such  high  lati¬ 
tudes  most  bodies  of  water  are  covered  with  Ice,  even  In 
summer . 

The  deviation  of  the  directions  determined  for 
July  13  Is  not  unreasonable,  considering  the  p»s.  Some  of 
the  microseisms  could  have  their  origin  on  the  eastern 
coast  of  U.  S.  A.  The  direction  to  the  southern  edge  of 
Hudson  Bay  is  a  =  +  8l° . 

The  directions  on  August  13,  6h  GMT  and  12h  GMT 
agrees  well  with  a  =  60°,  which  is  the  azimuth  of  the 
Hudson  Bay  low  and  the  Labrador  low.  At  9h  GMT  the  source 
seems  to  be  further  south.  The  determination  on  August  l4 
deviates  70°,  but  for  p  =  3  such  an  error  may  occur. 

The  directions  on  August  28  agree  well  with  the 
disturbed  weather  in  Greenland  and  in  the  Canadian  Arctic. 

4 . 5e  The  Distribution  of  the  Directions  of  Approach.  A 

Selective  Microseismic  Barrier  in  Western  Canada. 

A  significant  tendency  for  the  3-6  second 
microseisms  to  arrive  from  the  octants  N,  NE,  and  E  was 
demonstrated  in  section  4.4a. 

This  feature  may  be  explained  by  assuming  a 
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barrier  for  such  a  microseisms  between  the  Pacific  and 
Edmonton.  None  of  the  directions  of  approach  determined 
for  microseisms  with  periods  between  3  and  6  seconds  fell 
into  the  sector  between  a  =  -71°  and  a  =  -136°,  in  which 
the  distance  to  the  Pacific  is  shorter  than  to  Hudson  Bay. 
Since  all  the  observations  with  periods  between  5  and  6 
seconds  originate  from  only  two  storms,  the  conclusion  of 
a  microseismic  barrier  in  western  Canada  has  been  restricted 
to  the  3-5  second  microseisms,  but  it  is  possible  that  the 
relative  rareness  of  5  -  6  second  microseisms  is  caused  by 
the  same  effect. 

A  large  attenuation  of  3  -  5  second  microseisms 
of  Pacific  origin  also  agree  with  the  discussion  in  section 
4.5b,  and  in  section  4.5c  it  is  shown  that  an  8  second 
microseismic  storm  of  Pacific  origin  was  recorded,  hence  the 
barrier  in  western  Canada  can  be  described  as  selective. 

The  interval  of  azimuths  pointing  towards  the 
Rocky  Mountain  system  is  from  a  =  +145°  to  a  =  -26°.  This 
189°  interval  contains  2  out  of  the  24  determinations  with 
3-6  second  periods.  One  of  these  directions  was  from 
July  2.  Its  deviation  from  the  margin  a  =  +145°  was  less 
than  thaprobable  error.  The  other  direction  was  from  July 
19;  it  had  p  =  1  and  the  smallest  AR  determined,  hence  its 
reliability  is  extremely  low.  Omission  of  two  points  could 
change  its  direction  to  north. 
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Thus  the  observations  are  consistent  with  the 
microseismic  barrier  along  the  Rocky  Mountains,  but  this 
barrier  could  be  anywhere  between  Edmonton  and  the  Pacific 
coast . 

The  velocity  of  surface  waves  with  periods  3-8 
seconds  have  been  investigated  by  many  authors;  examples 
are  Gutenberg  (1955)*  Lehmann  and  Ewing  (i960)  and  Oliver 
(1962).  Velocities  were  generally  between  2  and  3-5  km/sec. 
Thus  wavelengths  for  4  second  microseisms  could  be  8  km, 
while  they  could  be  30  km  for  8  second  microseisms. 

The  selective  microseismic  barrier  found  could 

be  explained  if  the  severe  lithological  anomalies  in 

Western  Canada  penetrate  to  the  depth  comparable  to  8  km, 

but  not  to  30  km.  Conclusions  of  this  kind  must  be  delayed, 

however,  until  the  nature  of  the  continental  waveguides 

is  settled,  but  it  can  be  stated,  that  a  final  model  of 

the  crustal  structure  of  Western  Canada  must  be  able  to 

explain  this  selective  barrier,  unless  surface  topography  can 
account  for  it. 

This  study  is  the  first  to  demonstrate  a  barrier 
effect  in  Western  Canada,  but  similar  conclusions  have 
been  obtained  in  U.  S.  A.  as  shown  by  the  following  quota¬ 
tion  from  Oliver  and  Ewing  (1958b): 

"Carder  found  that  all  but  two  stations  in  the 

United  States  recorded  microseisms  of  4  to  7  seconds  from 
a  storm  off  the  east  coast.  The  exceptions  were  Berkeley 
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and  Ukiah .  Carder  attributed  this  anomaly  to  the  effect 
of  either  Sierra  Nevada  or  the  Central  Valley  of  California 

.  Carder  also  found  a  suggestion  of  better  micro- 

seismic  transmission  across  the  eastern  low  lands  than  across 
the  Rocky  Mountain  system." 

The  fact  that  no  5  -  9  second  microseisms  were 
found  to  arrive  from  the  east  is  considered  an  effect  of  the 
summer  season,  where  generation  of  such  microseisms  is  less 
prominent  than  in  winter.  The  continental  path  between 
Edmonton  and  the  east  coast  is  expected  to  transmit  such 
microseisms  efficiently;  at  least  the  8-9  second  periods, 
since  Donn  (1954)  found  Pacific  microseisms  with  such 
periods  at  the  Palisades  with  good  directional  characteristics. 

4.5f  Validity  of  the  Basic  Assumptions 

The  basic  assumptions  of  Jensen’s  method  are 
stated  in  section  2.5* 

The  first  assumption  requires  that  the  particle 
motion  of  the  Rayleigh- type  waves  is  retrograde.  No 
authors  have  reported  prograde  particle  motion  for  periods 
longer  than  3  seconds,  and  the  discussion  in  section  1.2 
shows  that  even  M2  waves  can  be  expected  to  have  retrograde 
particle  motion  for  periods  longer  than  3  seconds.  The  same 
discussion,  however ,  would  make  it  reasonable  to  expect 
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prograde  M2-waves  if  the  pertinent  crustal  model  has  an 
unusual  low  rigidity. 

Since  station  Edmonton  is!. located  on  unconsolidated 
shales,  on  a  thick  sedimentary  cover,  it  may  appear  necessary 
to  check  if  the  particle  motion  at  Edmonton  for  surface 
waves  with  periods  longer  than  3  seconds  is  retrograde. 

Such  a  check  is  provided  by  the  agreement  between 
the  observations,  and  the  cases  of  disturbed  weather  in 
Hudson  Bay,  the  Arctic  and  in  the  Pacific  as  described  in 
sections  4.5c  and  4.5d.  In  particular  the  observations 
from  July  6  with  p  =  9  and  p  =  10  occurred  with  disturbed 
weather  in  the  Arctic  and  quiet  conditions  at  the  Californian 
coast,  so  that  the  source  of  these  microseisms  can  be 
considered  well  established. 

No  research  was  conducted  specifically  to  test  the 
validity  of  the  second  basic  assumption,  but  the  good  agreement 
between  observations  with  the  large  p  value  and  the  weather 
situation  favors  the  adoption  of  this  assumption. 

Conclusions 

As  a  result  of  the  discussions  in  this  chapter 
it  can  be  concluded: 

1.  A  microseismic  barrier  in  western  Canada  attenuated 

the  3-5  second  microseisms  of  Pacific  origin,  so  that 
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it  was  never  observed  to  dominate  over  the  other 
microseisms,  but  a  storm  of  8  second  microseisms  of 
Pacific  origin  was  demonstrated. 

2.  The  meteorological  conditions  associated  with  the  8 
second  microseismic  storm  favors  the  conclusion  that  it 
is  caused  by  the  interference  of  incoming  and  reflected 
swell  at  the  Pacific  coast  in  agreement  with  Longuet- 
Higgins*  theory. 

3.  The  period  range  observed  for  microseismic  storms  of 
non-Pacific  origin  agrees  with  the  range  predicted 
from  the  wind  velocities  using  equation  (4-1)  and 
Longuet-Higgins 1  theory. 

4.  The  5-6  second  microseisms  observed  at  Edmonton 
during  July  and  August  occurred  in  two  storms  andf 
originated  in  Arctic  waters. 

5.  The  most  frequent  approach  for  3-5  second  micro¬ 
seisms  was  NE  and  E. 

6.  Good  correlation  was  found  between  the  directions  of 
approach  determined  for  microseismic  storms,  and 
the  directions  expected  from  the  disturbed  weather  at 
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V.  STUDY  OF  THREE  CYCLES  PER  SECOND  MICROSEISMS 

5 . 1  Introduction 

While  2  cps  microseisms  are  found  quite  persistently 
in  many  areas  of  the  United  States  (Wilson,  Vesiac  Staff, 

1962)  a  similar  spectral  peak  has  not  been  found  in  Alberta. 
During  the  1962  field  season  no  2  cps  wavelets  were  seen  on 
the  field  records,  but  on  the  University  of  Alberta  Observatory 
short  period  records  it  does  show  up  occasionally  (visual 
observation)  .  Frequencies  around  3  cps,  however,  show  up 
quite  consistently  (but  not  everywhere  at  all  times),  and 
it  is  concluded  that  in  Alberta  the  peak  of  the  power 
spectrum  tends  to  be  at  3  cps  rather  than  2  cps. 

In  analogy  with  chapter  IV  the  first  sections  of 
this  chapter  will  present  the  experimental  data  obtained 
at  University  of  Alberta,  while  the  comparison  with  work 
done  elsewhere,  the  interpretation  and  the  final  conclusions 
will  follow  in  later  sections. 

While  the  field  records  showed  that  some  locations 
in  general  had  a  much  larger  3  cps  amplitude  than  other 
locations,  the  following  section  describes  a  study  showing 
that  at  one  given  location  the  average  amplitude  is  much 
larger  at  some  times  than  at  other  times. 
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5 . 2  Comparison  Between  Observatory  Records  and  Local  Wind 

On  the  short-period  observatory  records  many 
individual  bursts  of  energy  of  a  few  minutes  duration  are 
present.,  some  of  them  having  a  very  regular  appearance  like 
3  cps  sinusoidal  waves.  Between  these  bursts  some  rather 
irregular  background  microseisms  are  always  present  and  have 
a  high  content  of  the  3  cps  frequency.,  though  the  2  cps 
frequency  is  also  present.  A  glance  at  a  few  records  reveals 
a  daily  variation  (indicating  more  quiet  conditions  at 
night),  and  this  was  originally  thought  mainly  to  be  due  to 
wind . 

Figure  8(a)  illustrates  this  relationship  by  com¬ 
paring  the  background  noise  with  the  wind  that  was  simultane¬ 
ously  observed  at  the  International  Airport.  The  amplitudes 
of  the  background  noise  were  measured  by  ignoring  the 
discrete  periods  of  high  activity,  caused  by  traffic  and 
farm  machinery.  The  amplitudes  were  measured  during  the 
period  June  12  -  30,  1963.  Surprisingly,  no  clear  correlation 
between  background  microseisms  and  wind  is  visible .  To 
investigate  whether  this  lack  of  correlation  could  be  due  to 
variation  of  wind  from  one  location  to  anothei}  figure  8(b) 
was  made.  The  University  of  Alberta  Observatory  is  12  miles 
from  the  International  Airport  and  the  Municipal  Airport  is 
17  miles  from  the  International  Airport.  A  clear  correlation 


' 


WIND  AT  MUNICIPAL  AIRPORT  (MPH)  TRACE  AMPLITUDE  (MM) 


76. 


WIND  AT  INTERNATIONAL  AIRPORT  {  MPH) 

FIG.  8.  June  12-30  1963,  one  dot  per  simultaneous  3* hourly 
observation  of : 

(a)  U.ofA.  observatory  SPZ  3cp$  background  noise 
versus  wind  velocity  at  International  Airport. 

(  b)  Wind  velocity  at  Municipal  Airport  versus  wind 
velocity  at  International  Airport. 
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is  demonstrated  between  the  simultaneous  winds  at  the  two 
airports,  and  although  the  scatter  is  appreciable,  the 
correlation  is  so  good  that  it  is  possible,  based  on  figure 
8(b),  to  make  a  reasonable  correction  to  the  instruments 
that  were  used:  The  points  in  figure  8(b)  do  not  group 
symmetrically  about  the  identity  line  (the  full  line  on 
figure),  as  they  would  if  the  frequency  of  occurrence  of  any 
given  wind  velocity  were  similar  at  the  two  places  and  if 
the  instruments  were  reading  correctly.  However,  the  distri¬ 
bution  is  symmetrical  about  the  dashed  line,  which  is  the 
identity  line  that  would  have  resulted  if  the  instruments 
at  the  Municipal  Airport  were  reading  as  little  as  1  mph 
more  and  the  instruments  at  the  International  Airport  were 
reading  1  mph  less  than  true  wind.  Such  small  differences 
may  be  difficult  to  detect  by  other  means. 

Now  assume:  1.  Differences  in  wind  velocities 
between  the  Edmonton  International  Airport  and  the  Municipal 
Airport  are  comparable  to  differences  in  wind  velocities 
between  the  International  Airport  and  the  University  of 
Alberta  Observatory. 

2.  Wind  is  the  major  source  of  the  2-3  cps 
background  noise  at  the  University  of  Alberta  Observatory. 

Then  the  spread  observed  in  figure  8(b)  is  too 
small  to  make  the  expected  correlation  between  wind  at  the 
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FIG.  9.  June  12-30  1963,  distribution  of  observations  as 
a  function  of  the  time  of  day: 

(a)  U.  of  A.  observatory  SPZ  3cps  background 
noi  se 

(  b)  Wind  velocities  at  International  Airport . 
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International  Airport  and  background  noise  at  the  observatory 
disappear .  Now,  since  no  significant  correlation  appears  in 
figure  8(a),  and  since  assumption  1.  is  reasonable  (the  dis¬ 
tances  between  the  two  pairs  of  locations  are  comparable)  it 
is  possible  to  conclude  that  2.  is  wrong. 

The  conclusion  is  substantiated  by  figure  9.  Figure 
9(a)  -shows  the  amount  by  which  the  distribution  of  this 
observed  2-3  cps  background  noise  varies  with  the  time  of 
day.  Figure  9(b)  shows  that  for  the  period  of  observation 
(June  12  -  30,  1963)  no  corresponding  variation  of  wind 
velocities  with  the  time  of  day  is  present.  This  lack  of 
variation  of  wind  between  day  and  night  is  surprising,  and 
may  be  special  for  the  period  considered,  but  it  proves  that 
local  wind  is  not  the  main  source.  Even  the  maximum  winds 
for  each  day,  the  velocity  and  the  hour  of  occurrence  of 
which  was  noted,  did  not  seem  to  have  had  any  effect  that 
was  visible  on  the  corresponding  seismograms. 

Thus  it  has  been  demonstrated  that  the  local 
winds  are  not  the  main  cause  of  the  2-3  cps  background 
microseisms  at  station  Edmonton.  Furthermore  a  clear 
daily  variation  of  these  microseisms  has  been  demonstrated. 
This  variation  will  be  discussed  further  in  section  5.6c. 
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5-3  Microseisms  from  Trans- Canada  Highway 
5 . 3a  Introduction 

This  section  will  show  that  highway  traffic  is  one 
of  the  sources  of  3  cps  microseisms.,  and  some  data  concerning 
the  attenuation,  phase  velocity  and  particle  motion  will  be 
presented . 

5 . 3^  The  Set-up 

Figure  10  shows  the  location  of  the  L- spread  and 
of  the  three- component  station  (Z-N-E).  Automobile  traffic 
did  not  produce  noticable  3  cps  microseisms.,  but  the  very 
large  trucks  (more  than  ten  tons)  caused  the  microseisms  to 
increase  by  a  factor  of  ten  or  more. 

Three  such  recordings  were  made  taking  note  of  the 
truck  positions  with  respect  to  the  numbered  telephone  poles 
shown  in  figure  10;  the  corresponding  trucks  will  be  called 
1  „  2,  and  3.  When  truck  1  was  passing,  only  four  vertical 
seismometers  were  in  operation  (  in  positions  corresponding 
to  trace  6,  7,  8,  and  9  on  figure  10) .  Truck  1  was  moving 
fifty-four  miles  per  hour,  while  trucks  2  and  3  were  moving 
about  forty  miles  per  hour. 
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O  TRACE  I 
O  TRACE  2 
O  TRACE  3 

O  TRACE  4 
TRACE  6 
TRACE  7 
TRACE  8 
•—TRACE  10,114.12  (  E ,  N  &  Z  ) 
TRACE  9 


I: 20000 
I  km 


Point  A  is  SWC.section 5,  Tp  18,  R  12,  W4. 


FIG.  10. 
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5 . 3c  The  Attenuation 

A  number  of  amplitudes  were  measured  in  order  to 
get  information  about  the  attenuation.  Amplitudes  from 
truck  1  were  multiplied  by  .27  to  make  them  comparable  to 
those  due  to  truck  2  and  truck  3.  This  factor  was  arrived 
at  by  plotting  the  amplitudes  on  a  logarithmic  scale  as  a 
function  of  distance  and  then  considering  thepattern  of  points 
produced  by  each  truck. 

In  figure  11  these  amplitudes  are  plotted  as  a 
function  of  distance  between  truck  and  seismometer.  Differ¬ 
ent  symbols  are  used  for  different  groups  of  traces  to  test 
whether  the  amplitudes  vary  smoothly  with  distance  or  if 
horizontal  variation  of  elastic  parameters  along  the  setup 
caused  amplitude  anomalies. 

The  points  in  the  range  650  -  750  meters  on 
figure  11  show  that  traces  8,  9,  and  12  consistently  have 
larger  amplitudes.  In  the  range  800  -  1100  meters  the 
amplitudes  on  trace  8,  9,  and  12  differs  so  much,  that  the 
points  encircled  on  the  figure  were  investigated.  It  then 
turned  out  that  these  points  correspond  exactly  to  all  the 
readings  on  traces  8,  9,  and  12  originating  from  truck 
positions  near  pole  0  and  pole  1. 

Two  different  explanations  can  be  offered  for  this 
effect:  First,  the  elastic  properties  of  the  ground  may 

possess  an  antinode  for  3  cps  vibrations,  making  excitation 
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FIG.  II  „  Ground  velocity  (v)  versus  distance  (r)  from  truck  on  Transcan  a  da  Highway. 
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of  these  vibrations  more  efficient.  Second,  the  effect  may 
simply  be  due  to  an  increased  noise  enhancement  associated 
with  the  slight  change  of  the  direction  of  the  truck  move¬ 
ment  (see  figure  10). 

Whatever  the  cause  for  the  increased  microseisms 
is,  it  seems  reasonable  to  reject  all  the  values  originating 
near  poles  0  and  1.  Only  truck  2  and  truck  3  were  recorded 
in  that  area,  and  hence  as  a  consequence  of  the  rejection 
of  the  points  it  was  found  that  the  factor  .27.,  used  to 
make  the  amplitudes  from  truck  1  comparable  to  those  of  the 
other  trucks,  had  to  be  changed  to  .17. 

By  doing  this  the  scatter  of  the  points  was  reduced 
considerably,  and  figure  12  shows  the  remaining  points  with 
the  ordinate  changed  from  v  to  v/r,  where  r  is  the  distance 
in  kilometers  and  v  the  ground  velocity  in  Angstroms  per 
second.  A  straight  line  on  the  diagram  will  correspond  to 
a  constant  value  of  the  attenuation  coefficient  a. 

Figure  12  shows  that  the  points  corresponding  to 
traces  1,  2,  3,  4,  6,  and  7  follow  the  same  attenuation 
pattern,  but  the  points  corresponding  to  the  amplitudes  on 
traces  8,  9,  and  12  are  scattered  at  a  significantly  higher 
level  (50^  -  100^  higher).  It  is  worth  noting  that  trace 
8,  9,  and  12  seismometers  were  all  within  150  meters  of  a 
shallow  pond  while  the  other  seismometers  were  on  the  flat 
and  dry  prairie  far  from  any  topographic  anomalies. 
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The  line  drawn  corresponds  to  a  =  3.8  km-1 

and  the  mean  frequency  Is  3.1  cps .  This  agrees  well  with 
Kisslinger  1 s  value,  4  to  6  km"1  at  3-5  to  5  cps,  found  at 
the  Suf field  Experimental  Station  (Jones,  1963) . 

The  points  corresponding  to  truck  positions  more 
than  two  kilometers  away  do  not  fit  the  attenulation  pattern, 
and  they  may  be  due  to  microseisms  from  other  sources. 

5 . 3d  The  Phase  Velocity 

In  order  to  determine  the  phase  velocity  of  the 
waves  and  investigate  the  spatial  coherency  of  the  waves, 
the  velocity  of  wave  crests  in  well  developed  groups  recorded 
from  truck  1  was  determined.  Similar  studies  using  truck  2 
or  truck  3  did  not  produce  meaningful  results,  because  the 
waves  were  more  irregular. 

The  average  value  for  the  difference  At  of  arrival 
times  on  trace  9  and  trace  6  was  determined  for  selected 
wavegroups.  The  truck  position  at  the  time  of  the  wave 
generation  was  calculated  and  the  difference  Ar  between 
distances  to  trace  6  and  trace  9  seismometers  was  determined. 
Positive  and  negative  Ar 1 s  occurred. 

The  fact  that  At  was  decreasing  smoothly  to  zero 
for  Ar  decreasing  to  zero  is  taken  as  a  sign  that  phase 
identification  is  correct. 
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The  mean  period  derived  from  59  measurements  is 
0.312  -  0.002  seconds ,  and  an  estimate  obtained  for  the 
spectral  width  parameter  6  is  0.09. 

The  narrowness  of  the  spectrum  makes  it  reasonable 
to  assume  that  the  apparent  wave  velocity  is  a  good  estimate 
of  the  phase  velocity. 

Assuming  a  rectilinear  wave  propagation  from  the 
source  to  the  receiver,,  Ar/At  will  be  an  estimate  of  this 
phase  velocity. 

The  estimate  of  the  phase  velocity  obtained  for 
truck  positions  between  pole  4  and  7  and  also  for  truck 
positions  between  pole  0  and  1  is  489  -  6  m/sec.  For  truck 
positions  between  pole  1  and  3*  however,  the  estimate  obtained 
is  350  -  50  m/sec. 

The  latter  estimate  of  the  standard  deviation  was 
calculated  from  only  four  observed  deviations  and  hence  it 
can  be  considered  fairly  uncertain.  The  difference  between 
estimates  of  phase  velocities  could  be  due  to  the  decrease 
of  accuracy,  when  the  At 1 s  and  Ar 1 s  are  small,  but  at  least 
an  apparent  shift  of  phase  velocity  has  been  found. 

5 . 3e  The  Particle  Motion 

The  particle  motion  was  studied  using  the  modified 
Jensen  method.  Due  to  the  low  sensitivity  of  the  Willmore 
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Mark  I  seismometer  the  method  could  only  be  applied  to  a 
limited  number  of  cases.  The  diagrams  suggested  either 
predominantly  prograde  or  a  mixture  of  equal  amounts  of 
prograde  and  retrograde  Rayleigh-type  motion. 

The  trajectory  of  the  particle  motion  for  a  case 
when  the  NE  diagram  did  not  show  any  significant  trend  was 
followed  for  three  consecutive  cycles.  The  particle  motion 

was  the  same  as  would  result  from  a  refracted  SV  wave?  The 

0 

radial  amplitude  was  85  A/sec,  the  transverse  amplitude  five 

o 

times  less,  and  the  vertical  amplitude  was  93  A/sec. 

5 . 4  Microseisms  from  Highway  2 

5 . 4a  Introduction.,  Comparison  with  Trains 

In  order  to  collect  more  data  on  the  particle 
motion  in  3  cps  surface  waves *  some  recordings  were  made 
near  highway  2.  The  site  selected  was  as  close  to  the 
Edmonton  seismic  observatory  site,  as  it  was  possible  to 
obtain  known  sources  of  3  cps  microseisms.  These  sources 
were  trains  and  trucks .  Comparison  of  the  amplitudes 
suggests  that  the  3  cps  energy  produced  by  a  heavy  truck 
compares  with  the  3  cps  energy  produced  by  one  car  of  a 
freight  train^  but  from  one  car  ol  a  passenger  train  the 
energy  is  less.  The  following  analysis  is  based  on  record¬ 
ings  of  trucks . 

i.e.  a  vertically  polarized  S  wave. 
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5  •  4  b  The  Set-up,  Description  of  the  Recordings 

The  setup  is  shown  in  figure  13.  The  first  three 
traces  form  a  tripartite  recording  with  a  distance  of  150 
ft.  between  the  seismometers.  At  the  same  location  the  N 
and  E  components  are  recorded. 

Figure  14  shows  a  section  of  one  of  the  recordings, 
when  the  corresponding  truck  was  at  position  A  on  figure  13. 
The  sinusoidal  appearance  of  the  oscillations  is  apparent, 
especially  on  traces  4  to  11. 

The  attenuation  of  the  amplitude  with  distance  is 
clearly  visible  on  figure  14. 

5.4  c  The  Phase  Velocity 

Using  the  same  method  as  in  section  5.3d,  the 
phase  velocity  was  determined  from  the  differences  of  arrival 
times  at  trace  4  and  trace  6  seismometers.  The  azimuth  of 
the  source  varied  from  +  94°  to  +  50°  and  the  results  were 
consistent.  Their  mean  value  is  j60  +  19  m/sec. 

On  the  section  shown  in  figure  l4  no  phase  veloci¬ 
ties  can  be  determined  from  wave  correlation  from  the  tri¬ 
partite  set-up  on  trace  1-3;  but  when  the  trucks  were  in  the 
azimuth  range  +  28°  to  +  37°;  the  appearance  of  these  traces 
was  quite  regular  and  the  vectorial  wavenumbers  were  deter- 
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FIG.  13.  Set-up  near  highway  2. 
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mined  from  the  apparent  velocities  in  the  directions  S  and 
W.  The  vectorial  wavenumbers  are  plotted  on  figure  15 . 

The  average  phase  velocity  determined  from  these 
wavenumbers  is  880  -  80  m/sec,  but  the  pattern  exhibited 
by  the  end  points  of  the  vectorial  wavenumbers  is  not  as 
expected,  the  pattern  shows  a  large  change  of  the  apparent 
phase  velocity  with  azimuth. 

5 . 4d  The  Particle  Motion 

The  particle  motion  was  determined  by  the  modified 
Jensen  method  and  the  results  are  shown  in  figure  16. 

Ref.  1  shows  a  very  clear  case  of  prograde  Rayleigh 
waves,  the  sector  in  which  the  truck  moved  during  the  time 
used  for  the  measurement  is  opposite  to  the  sector  containing 
most  points.  The  distance  to  the  source  is  about  400  meters. 
The  transverse  motion  is  small. 

Ref.  2  is  derived  from  another  truck  having  a 
similar  range  of  positions  as  the  one  shown  on  Ref.  1.  The 
conclusion  is  the  same:  Prograde  Rayleigh  waves  with  a 
small  transverse  component.  Ref.  1  and  2  appear  consistent 
with  a  10  -  20  degree  lateral  refraction  of  the  waves. 

Ref.  3  shows  dominant  retrograde  Rayleigh  waves. 

The  distance  to  the  source  is  about  500  meters. 

Ref.  4  shows  dominant  retrograde  Rayleigh  waves 
from  a  source  700  meters  away. 
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FIG*  15  .  Vectorial  wavenumbers,  determined  from 
the  tripartite  set-up  shown  on  fig.  13. 
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NE  diagrams  used  "to  determine  the  particle  motion  of  3  cps 
microseisms  from  Highway  2.  (See  page  13  and  page  92) . 
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5 • 5  Regional  Microseism  Anomalies  in  Alberta 

The  analysis  in  section  5-3c  demonstrated  a 
regional  anomaly  of  the  amplitudes  of  3  cps  microseisms. 
Figure  17  shows  a  much  larger  regional  anomaly.  Traces 
7,  8,  and  9  have  a  maximum  amplitude  of  7  mm,  an  attenuation 
setting  of  6,  and  a  high  cut  filter  setting  of  8,  hence 

o 

from  (3-15)  the  ground  velocity  is  350  A/sec.  On  most  other 
traces  the  amplitude  is  1  mm  or  less  which  corresponds  to 

o 

50  A/sec.  This  anomaly  cannot  be  explained  from  the  topo¬ 
graphy.,  as  the  corresponding  seismometers  were  placed  in 
flat  irrigated  farm  land.* 

The  spread  was  on  the  north  line  of  section  33^ 

Tp  17,  R  l4,  W  4  ore  mile  from  Lake  Newell.  The  trace  one 
seismometer  was  on  the  west  end.  The  wind  speed  was  l6 
knots,  gusting  to  25  knots.  The  appearance  of  the  individual 
traces  continued  to  be  the  same  for  all  recordings  made. 

One  recording  at  the  shore  of  McGregor  Lake 
yielded  3  cps  microseisms  with  a  ground  velocity  of  950 
A/sec.  Two  foot  waves  with  a  period  of  2  seconds  were 
observed  in  the  lake,  but  shorter  wavelengths  were  super¬ 
imposed  on  these. 

Another  region,  where  unexplained  anomalously 
large  3  cps  amplitudes  were  observed,  was  one  mile  south  of 
Brant,  Tp  l8,  R  26,  W  4.  Here  the  winds  were  calm  and  the 


The  14  cps  microseisms  on  traces  2,  3;  4,  5  is  caused  by 
an  irrigation  canal  perpendicular  to  the  spread. 
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amplitude  variation  with  time  was  noted.  The  amplitude  of 
the  microseisms  varied  over  a  wide  range  in  a  few  minutes 
time;  this  variation  could  correlate  with  distant  highway 
and  railroad  traffic;  but  insufficient  information  is  avail¬ 
able  . 

5 . 6  Interpretation  and  Discussion 

It  is  not  possible  to  give  a  unique  interpretation 
of  all  the  data  presented  in  this  chapter,  because  often  a 
number  of  different  mechanisms  could  explain  the  same  data. 

5 . 6a  Attenuation,  Anomalous  Propagation 

The  attenuation  coefficient  found  near  the  Trans- 
Canada  Highway  is  3.8  km-1  while  a  tentative  value  found 
near  Highway  2  is  0.8  km”1;  in  both  cases  the  amplitudes 
would  become  insignificant  after  a  few  miles,  and  this 
agrees  with  the  observations .  It  also  explains  the  fact 
that  no  train  microseisms  could  be  identilied  on  the  short 
period  seismograms  obtained  from  station  Edmonton,  7  miles 
from  the  railroad  track.  Assuming  that  the  group  velocities 
in  the  areas  with  the  values  0=3. 8km  and  0=0. 8km  are 

0,3  km/sec  and  0.5  km/sec  respectively,  the  corresponding 

Q  values  will  be  8  and  24. 

These  values  seem  far  too  low  to  represent  energy 
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dissipation,  but  they  could  in  pant  be  due  to  loss  of  elastic 
energy  to  more  deep  seated  layers. 

Frantti  (1963)  suggested  that  a  substantial  portion 
of  the  ambient  high  frequency  earth  noise  originates  at  the 
earth* s  surface  and  follows  a  near  surface  propagation  path. 

Combining  Frantti* s  suggestion  with  the  Q,  values 
observed  above  it  is  suggested  that  3  cps  seismic  energy 
is  transferred  from  the  surface  to  a  near  surface  waveguide. 

That  this  waveguide  is  near  the  surface  was  con¬ 
firmed  by  B.  Isacks  (personal  communication),  who  reported 
that  in  the  recordings  made  by  the  Lamont  Geological  Obser¬ 
vatory  in  a  mine,  the  1  cps  component  had  similar  amplitudes 
at  the  surface  and  in  the  mine  at  2000 feet,  while  the  3  cps 
component  was  smaller  in  the  mine.  E.  J.  Douze  (personal 
communication)  using  a  deep  well  seismometer,  found  the 
3  cps  peak  to  disappear  with  depth. 

The  regional  anomalies  of  3  cps  microseisms 
reported  in  section  5*5  could  be  due  to  outcroppings 
of  this  waveguide  (i.e.  antinodes  at  the  surface). 

An  interesting  example  of  anomalous  behavior  of 
3  cps  wave  propagation  was  given  by  R.  A.  Page  (personal 
communication),  who  kindly  made  copies  of  such  seismograms 
available  for  me: 

Whenever  a  large  freight  train  was  passing  the 
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seismic  observatory  at  Sterling  Forest,  New  York,  two  discrete 
pulses  of  3  cps  microseisms  were  observed.  Duration  of  the 
first  pulse  was  about  four  minutes  with  a  maximum  ground 

o 

velocity  of  4000  A/sec,  then  five  minutes  of  usual  background 

o 

microseisms  with  amplitudes  less  than  1000  A/sec  and  then  a 
two  minute  pulse.  The  closest  approach  of  the  railroad  track 
from  the  vault  was  about  6.5  km.  The  phenomenon  is  currently 
being  studied  by  B.  I  sacks  and  R.  Page  at  Lam on t . 

By  comparison  with  the  University  of  Alberta 
findings  it  is  the  two  peaks,  rather  than  the  quiet  period 
in  between,  that  can  be  considered  anomalous.  These  peaks 
could  be  explained  by  two  outcroppings  of  the  waveguide 
along  the  railroad  track,  and  a  third  one  at  the  observatory. 

5 . 6b  Nature  of  the  Waves 

The  microseisms  discussed  in  section  5.3  and  5.4 
produced  by 

were  notAair- coupled  waves,  because  the  phase  velocity  was 
considerably  larger  than  the  speed  of  sound. 

The  particle  motion  was  most  often  a  clear  case 
of  either  retrograde  Rayleigh  type  waves,  most  likely  M-^, 
or  prograde  waves,  most  likely  Mg.  In  a  few  cases  the  NE 
diagrams  could  be  explained  by  a  mixture  of  prograde  and 
retrograde  Rayleigh  waves,  having  similar  amplitudes. 
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In  one  case  the  particle  motion  could  be  explained 
by  assuming  refracted  SV  waves.,  but  the  same  particle  motion 
can  be  explained  as  Interference  between  prograde  and  retro¬ 
grade  waves.  Indeed  any  Lissajou  figure  can  be  explained 
this  way:  Choose  a  line  in  the  plane  of  the  particle  motion. 
Assume  that  the  projection  of  the  prograde  and  retrograde 
component  on  this  line  are  of  equal  amplitude  and  l80°  out 
of  phase.  Then  the  particle  motion  will  be  rectilinear  on 
a  line  perpendicular  to  the  first.  Hence  by  superposition 
of  two  waves  the  particle  motion  can  degenerate  to  a  line  of 
any  orientiation,  and  for  other  amplitude  ratios  it  may 
become  an  ellipse  with  any  inclination  of  the  axis. 

5 . 6c  Sources  of  3  cps  Microseisms 

Freight  trains  and  heavy  trucks  were  demonstrated 
to  be  sources  of  3  cps  microseisms.  Another  source  of 
3  cps  microseisms  was  reported  by  Y.  Rochard  (Vesiac  Staff, 
1962,  chapter  15)  who  found  slightly  perceptible  3  -4  cps 
microseisms  from  a  stone  quarry  2  km  away . 

Farm  machinery  produced  a  high  level  of  3  cps 
microseisms  at  the  University  of  Alberta  observatory,  when 
such  machinery  was  operating  in  the  adjacent  tields. 

The  measurements  made  near  lakes  disturbed  by 
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winds  suggest  that  some  3  cps  microseisms  could  be  produced 
here,  but  evidence  is  Inconclusive. 

An  unsettled  question  is  the  nature  of  the  source  of 
the  2-3  cps  background  microseisms.  As  demonstrated  in 
section  5*2  they  show  a  pronounced  daily  variation.  The 
attenuation  pattern  observed  on  the  surface  suggests  that  the 
decrease  in  human  activities  at  night  cannot  account  for 
this  daily  variation,  but  a  stronger  suggestion  of  non- 
cultural  sources  was  obtained  from  P.  L.  Gouin  (personal 
communication) : 

"The  seismic  observatory  recordings  at  Addis 
Abeba  had  a  high  level  of  1.5  -  10  cps  microseisms,  which 
invariably  decreased  at  night  by  a  factor  of  5  or  so.  To 
filter  out  these  disturbing  microseisms  the  1/4  sec. 
galvanometers  were  replaced  by  2  sec.  galvanometers,  but 
the  effect  was  still  visible  with  a  smaller  trace  amplitude. 
These  microseisms  could  not  be  attributed  to  man-made 
sources;  because  the  heavy  traffic  took  place  at  night, 
where  the  microseisms  were  at  a  minimum.  Neither  could  they 
be  attributed  to  general  background  noise  adding  up  from  a 
large  number  of  man-made  sources  and  travelling  through 
some  waveguide ;  this  is  ruled  out  because  in  Ethiopia  no 
industry  is  present  and  the  traffic  is  limited  to  a  few 
lines  radiating  from  Addis  Abeba. 
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The  daytime  Increase  occurred  whether  the  sun 
was  shining  or  not.,  and  this  Increase  took  place  within 
half  an  hour  around  sunrise  and  died  out  within  half  an 
hour  around  sunset.  No  striking  correlation  with  the  ampli¬ 
tude  of  the  dally  variation  on  the  magnetograms  were  noticed” . 

This  information  gives  additional  support  to 
Gherzi’s  inference  from  the  Brebeuf  seismograms,  that  some 
short  period  noise  is  not  of  cultural  origin.  Gouin’s 
finding  makes  it  doubtful  that  the  daylight  heating  is  the 
source  but  still  makes  it  reasonable  indirectly  to  suspect 
the  sun  as  the  source . 

D.  E.  Smylie  (personal  communication)  suggested 
tentatively  that  such  microseisms  might  be  due  to  vibrations, 
set  up  in  the  ionosphere  by  solar  winds  or  by  other  electro¬ 
magnetic  disturbances , which  then  propagate  with  acoustic  speed 
to  the  ground . 

Another  suggestion  that  some  3  cps  microseisms 
are  of  atmospheric  origin  was  given  by  Walsh  (1955)*  He 
studied  such  microseisms  near  Saint  Louis  and  found  that 
they  occur  much  more  often  by  day  than  by  night  and  that 
they  correlate  well  with  microbarographic  activities . 

Conclusions 

1.  Heavy  trucks  and  freight  trains  are  sources  of  3  cps 
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microseisms;  at  distances  more  than  800  meters  the 
3  cps  peak  is  not  disturbed  by  higher  frequencies, 
but  at  closer  range  the  higher  frequencies  occur. 

2.  In  most  cases  these  microseisms  are  damped  out  within 

a  few  miles,  but  in  special  cases  efficient  propagation 
for  longer  distances  occurs. 

3.  These  3  cps  microseisms  are  not  aircoupled  waves. 

4.  They  contain  prograde  and  retrograde  Rayleigh  type  waves. 

5.  In  two  cases  an  apparent  change  of  phase  velocity  with 
azimuth  was  observed. 

6.  Industry  and  farm  machinery  are  sources  of  3  cps  micro- 
seisms  . 

7.  The  2-3  cps  background  microseisms  at  station  Edmonton 
have  a  clear  daily  variation. 

8.  This  variation  is  not  necessarily  due  to  human  activities 
it  is  not  due  to  wind,  but  could  be  of  atmospheric  origin 
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VI.  FOUR  TO  FIFTY  CYCLES  PER  SECOND  MICROSEISMS 

6 . 1  Introduction 

A  detailed  study  of  the  sources  of  these  micro¬ 
seisms  was  not  undertaken  but  some  obvious  sources  were 
located:  Machinery,  traffic  and  wind. 

6.2  Amplitudes 

The  smallest  amplitudes  encountered  by  University 
of  Alberta  were  recorded  near  Lake  Superior  in  a  tree- covered 
area  with  the  seismometers  on  bedrock. 

Recordings  of  microseisms  from  this  location  were 
kindly  supplied  to  me  by  D.  Robertson.  Some  oscillations 
correlated  from  trace  to  trace  and  could  be  identified  as 
electrical  noise.  Others  could  be  either  due  to  electrical 
noise  or  to  microseisms.  By  treating  them  as  microseisms 
it  was  possible  to  arrive  at  some  maximum  values  for  the 
ground  velocities.  The  actual  velocities  may  be  even  smaller. 

O 

Typical  values  were:  20  A/sec  at  5  -  10  cps, 

30  A/sec  at  20  cps.  These  values  are  slightly  larger  than 
the  lowest  values  reported  by  Frantti,  Willis  and  Wilson 
(1962)  in  their  survey  of  the  amplitudes  of  microseisms. 
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6.3  Particle  Motion 

Some  recordings  were  made  at  University  of  Alberta 
campus  to  study  the  particle  motion  in  microseisms  generated 
by  machinery. 

At  night  time  it  was  possible  to  record  unperturbed 
7  1/2  cps  microseisms  generated  by  the  University  Power  Plant. 
The  particle  motion  was  determined  by  the  modified  Jensen 
method.  The  resulting  NE  diagram  is  shown  in  figure  18(a) 
and  the  wave  numbers  determined  from  a  tripartite  setup  at 
the  same  location  is  shown  on  figure  18(b) . 

Some  lateral  refraction  of  the  waves  is  apparent,, 
and  their  particle  motion  is  prograde  with  a  small  trans¬ 
verse  component .  This  transverse  component  seems  to  be  out 
of  phase  but  coupled  to  the  vertical  component,  so  that  the 
same  direction  of  transverse  velocity  is  present  whenever  Z 
is  at  a  maximum.  This  coupling  may  be  due  to  the  fixed 
position  of  the  source  and  the  receiver  and  to  the  periodic 
nature  of  the  source.  The  waves  cannot  be  called  pseudo- 
Rayleigh  waves,  because  the  transverse  and  the  vertical 

components  are  out  of  phase. 

Figure  18(c)  is  an  NE  diagram  from  a  recording  of 
12  cps  microseisms  generated  by  the  air  conditioning  machine 
in  the  Physics  Building.  The  particle  motion  is  predominantly 

o 

retrograde  and  the  particle  velocity  3800  A/sec. 


This 
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NE  diagrams  used -to  determine 
FIG. 18.  the  particle  motion  of  7  cps 
and  12  cps  microseisms . 


a)  ;  prograde 
e)  :  retrograde 


See  page  13  and  page  105. 
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microseism  was  recorded  close  to  the  source,  while  the 
7  l/2  cps  was  recorded  500  meters  from  the  source.  In 
both  cases  the  particle  motion  agrees  with  Rayleigh- type 
surface  waves,  but  only  in  the  second  case  is  the  funda¬ 
mental  mode,  M-p  possible.  In  the  first  case  the  prograde 

►v,  e  t  t  o  H. 

particle  is  due  to  M2  waves  (Sezawa  waves)  or  still  higher 


mode  surface  waves 
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VII.  SUGGESTIONS  FOR  FUTURE  RESEARCH 

From  the  limitations  quoted  in  the  conclusions 
of  chapter  IV  to  VI  many  possibilities  for  further  research 
are  obvious;  for  instance  more  statistics  about  the  directions 
of  approach  at  Edmonton  and  about  their  dependence  on  the 
season  would  be  of  interest. 

Furthermore  from  the  discussion  in  section  1.2 

how 

it  appears  desirable  to  investigate  pseudo-Rayleigh  waves 
occur.  This  kind  of  wave  is  a  special  case  of  a  coupling 
between  Rayleigh  and  Love  waves,,  namely  the  case  of  0°  or 
180°  phase  difference  between  the  vertical  and  the  trans¬ 
verse  particle  motion.  A  coupling  between  second  mode 
Rayleigh  and  Love  waves  can  be  expected  from  the  similarity 
of  the  dispersion  curves  (Oliver  and  Ewing,  1958a),  but  only 
in  special  cases  would  pseudo-Rayleigh  waves  occur  by  this 
mechanism . 

If  pseudo-Rayleigh  waves  or  coupled  waves  with  a 
constant  phase  difference  not  equal  to  90°  occur,  the 
particle  velocity  vector  at  times  where  Z  has  a  zero-up- 
crossing  will  have  a  constant  direction  perpendicular  to  the 
direction  of  approach. 

Consequently  in  cases  of  large  p  values  it  is 
s ted  to  read  the  slopes  at  times  where  Z  has  a  zero  up 
crossing  and  plot  -^r)  i-n  an  NE  diagram.  If  these 
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points,  when  sufficiently  separated  in  time,  fall  in  one 
direction,  the  Rayleigh  and  Love  waves  must  he  coupled. 

If  then  the  angle  between  the  line  through  the  points  and 
the  direction  of  approach  is  significantly  different  from 
90°,  the  waves  would  not  be  pseudo-Rayleigh  waves  but  an¬ 
other  kind  of  coupled  waves.  In  such  cases  the  perpendicular 
to  the  transverse  motion  in  the  half  plane  containing  the 
previously  determined  direction  of  approach  would  be  an 
improved  estimate  of  the  direction  of  approach,  because  if 
Z  =  0  and  the  motion  would  have  no  longitudinal  component. 

The  improvement  of  the  estimate  is  only  valid,  however,  if 
the  vertical  component  is  caused  by  one  kind  of  wave  coming 
from  one  direction. 

If  the  points  plotted  fall  evenly  in  two  directions 
l80°  apart,  the  transverse  component  must  be  due  to  uncoupled 
Love  waves,  and  the  best  line  through  the  points  may  serve 
the  purpose  of  getting  another  estimate  of  the  direction  of 
approach;  the  indeterminacy  of  l80°  is  not  important  once 
the  approximate  direction  is  known. 


r 
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APPENDIX  A 

RAYLEIGH  WAVES  IN  A  SOLID  LAYER  OVER  A  SOLID  HALF  SPACE 

A . 1  General  Discussion 

While  the  theory  of  Rayleigh  waves  in  a  homogeneous, 
isotropic  and  solid  half  space  is  rather  simple,  there  is  a 
considerable  increase  in  algebraic  difficulties  whenever  an 
additional  horizontal  layer  is  added  to  the  model.  Consequently 
most  calculations  are  performed  on  a  two  layer  case,  i.e.  a 
model  consisting  of  a  homogeneous,  isotropic  and  solid  half 
space  in  welded  contact  with  a  homogeneous,  isotropic  and 
solid  layer  of  constant  thickness.  In  this  appendix  only 
two  layer  cases  will  be  considered;  the  primed  quantities 
will  refer  to  the  superficial  layer  and  the  unprimed  quantities 
will  refer  to  the  half  space. 

The  following  symbols  will  be  used: 

L  =  wavelength 
k  =  wave  number  =  2tt/L 
H  =  layer  thickness 

v  =  S-wave  velocity 
s 

=  Lame's  elastic  parameters,  \x  =  rigidity 
p  =  density 

C  =  Rayleigh  wave  velocity  in  a  half  space 
R 

" » "  refers  to  upper  medium  . 
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By  solving  the  period  equation  (e.g.  Ewing,  Jardetzky 
and  Press,  1957  ^  p.  191)  derived  from  the  boundary  conditions 
on  the  surface  and  the  interface  it  is  possible  to  calculate 
the  phase  velocity  as  a  function  of  the  dimensionless  para¬ 
meter  kH.  Several  solutions  are  possible,  corresponding  to 
different  modes.  Among  these  only  the  fundamental  Rayleigh 

mode  will  have  a  phase  velocity  that  approaches  CL 1  when 

n 

kH  ^  oo.  The  phase  velocities  in  the  other  modes  will  in 

general  approach  v  ’  for  kH  -*■  <». 

s 

If  the  rather  stringent  condition  for  the  existence 

of  Stoneley  waves  is  satisfied,  there  will  be  an  additional 

mode  for  which  the  phase  velocity  approaches  the  speed  of 

Stoneley  waves  as  kH  ->  oo.  The  condition  for  the  existence 

of  Stoneley  waves  is  shown  by  Ewing,  Jardetzky  and  Press 

(1957,  p.  113)  for  two  different  cases;  the  condition  is  that 

the  difference  between  the  S  wave  velocities  v  and  v  1  is 

s  s 

fairly  small. 

An  infinity  of  additional  modes  exists  with  phase 
velocities  approaching  vg  1  when  kH  ->00.  These  are  the  higher 
shear  modes,  each  of  which  has  a  cutoff  frequency,  correspond¬ 
ing  to  a  least  value  of  kH.  Below  this  cutoff  frequency 
unattenuated  propagation  does  not  occur. 

The  higher  shear  modes  are  numbered  so  that  the 
first  shear  mode  in  general  has  the  lowest  cutoff  frequency 
and  the  lowest  phase  velocities.  At  the  cutoff  frequency 
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the  phase  velocity  is  equal  to  v  . 

s 

In  the  model  computed  by  Tolstoy  and  Usdin  (1953, 
p.  860)  tne  first  three  shear  modes  are  computed  by  the 
method  of  ray  superposition;  in  their  notation  the  first 
three  shear  modes  are  M^,  M.^,  and  .  The  distinction 
between  the  modes  Mln  and  M2n  originates  from  the  theory  of 
vibrations  in  elastic  plates,  where  corresponds  to  the 
symmetric  and  to  the  antisymmetric  vibrations.  According 

to  J.  Oliver  and  L.  R.  Sykes  (personal  communication)  this 
distinction  is  not  convenient  for  waves  in  a  layered  half 
space,  and  the  waves  are  far  from  being  either  symmetric  or 
antisymmetric . 

In  the  model  used  by  Tolstoy  and  Usdin,  the 
is  the  first  shear  mode;  this  mode  is  called  by  many 
authors,  including  Oliver  and  Ewing  (1957).  The  question 
of  the  nomenclature  of  the  higher  mode  Rayleigh  waves  is 
not  settled,  however;  sometimes  the  numbering  of  the  shear 
modes  is  arbitary  when  phase  velocity  curves  intersect.  A 
classification  based  on  the  number  of  nodal  planes  for  the 
horizontal  or  vertical  motion  is  not  immediately  recommend - 
able,  because  no  general  simple  rule  is  known  that  gives  the 
number  of  such  nodal  planes  in  a  given  shear  .mode.  In  the 
fundamental  Rayleigh  mode  there  is  one  nodal  plane  for  hori¬ 
zontal  motion  and  no  nodal  plane  for  the  vertical  motion. 
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A.  2  Particle  Motion  and  Rayleigh's  Constant  as  an  Exploration 


Tool 


As  mentioned  in  section  1.2  prograde  particle 
.motion  on  the  surface  will,,  in  the  mode.,  usually  be  associ¬ 
ated  with  a  large  ratio  ,  while  smaller  ratios  will  be 

associated  with  retrograde  motion.  When  prograde  particle 
motion  is  observed,  it  can  be  concluded  that  the  mode  observed 
is  not  the  fundamental  Rayleigh  mode.  Assuming  the  mode  to 
be  it  can  then  be  concluded  that  the  ratio  pi/pt 1  is  rather 

large.  Further  conclusions  can  be  reached  by  considering  the 
Rayleigh's  constant,  i.e.  the  horizontal  amplitude  divided 
by  the  vertical  amplitude.  For  convenience  we  will  define 
the  Rayleigh  constant  to  be  positive  for  retrograde  and 
negative  for  prograde  particle  motion. 

Lee  (1932  and  1934)  suggested  that  Rayleigh's 
constant  provides  a  useful  tool  to  get  information  about 
the  superficial  layers  of  softer  sediments.  Since  this 
information  can  be  deduced  without  knowing  the  distance  to 
the  source,  the  microseisms  may  provide  the  Rayleigh  waves 
needed . 

The  Rayleigh  constants  for  M1  and  M 2  waves  will 
both  be  functions  of  L/H,  |a/p',  p/p'.  In  general  the 
variation  of  the  rigidity  contrast  is  much  larger  than  the 
variation  of  the  density  contrast,  and  in  the  models  chosen 
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by  different  authors  the  density  contrast  p/p'  has  always 
been  close  to  1  while  the  rigidity  contrast  \i/\il  has  varied 
between  2  and  oo.  Table  V  shows  some  values  of  Rayleigh's 
constant  for  the  M^mode  while  Table  VI  shows  some  values  of 
this  constant  for  the  M2  mode.  Not  all  the  values  given 
originate  from  just  one  source:  3  points  of  Sezawa  and 
Kanai '  s  (1940)  calculations  agreed  well  with  the  curves 
given  by  Suzuki  (1933),  and  Oliver  and  Ewing  (1957)  checked 
3  of  their  results  against  Kanai 's  earlier  values,  calculated 
with  the  same  parameters. 

The  usefulness  of  Rayleigh's  constant  as  a  tool  to 
obtain  information  about  parameters  of  crustal  layers  was 
questioned  by  Banerji,  who  wrote  that  Lee's  analysis  for 
layers  that  are  of  extent  small  in  comparison  with  the 
wavelength  of  the  Rayleigh  waves  appear  to  be  unsound.  The 
validity  of  Banerji's  objection  will  depend  on  the  layer 
thickness  and  on  the  rigidity  contrast.  As  Table  V  shows, 
a  rigidity  contrast  of  1:5  and  a  wavelength  20  times  the  layer 
thickness  will  give  a  Rayleigh  constant  of  0.99,  which  is 
appreciably  higher  than  the  value  0.68  approached  for  ^  -*  00. 

Since  the  Rayleigh  constant  depends  on  many 
parameters,  and  since  a  two  layer  model  is  not  always 
detailed  enough  to  explain  experimental  data,  the  Rayleigh 
constant  alone  will  provide  very  incomplete  information 
about  crustal  structure;  it  may  provide  the  check,  however. 
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TABLE  V 


Rayleigh's  Constant  in  the  W  Mode  with 
Poisson's  condition  satisfied  in  both  media. 


Source  of 

p/p ' 

i  M-/U  ' 

... 

1 

Rayleigh's 

information 

L/H 

constant 

Suzuki  (1933) 

1 

2 

1 

0.67 

1 

2 

2 

0.76 

1 

2 

15 

0.8l 

1 

2 

20 

0.78 

1 

4 

1 

0.67 

1 

4 

2 

0.70 

1 

4 

9 

1 . 14 

1 

4 

20 

0.96 

1 

5 

1 

0.67 

1 

5 

3 

0.48 

1 

5 

6 

1.19 

- 

1 

5 

8 

1.30 

1 

5 

20 

0.99 

Sezawa  and 

1 

20 

2 

0.64 

Kanai  ( 1940) 

1 

20 

5 

2.00 

1 

20 

16 

2.63 

1 

00 

2 

0.59 

1 

00 

6 

0.34 

1 

00 

9 

0 

no 

uo 

Model:  '  t 

H  p»  \i '  =  A* 

it 


P 


A 
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TABLE  VI 


Rayleigh's  Constant  in  the  M^  Mode  with 
Poisson's  Condition  Satisfied  in  both  Media. 


Source  of 

information 

p/p' 

( 

I 

ia/U  ' 

L/H 

Rayleigh ' s 

Constant 

Oliver  and 

i — 1 

i — 1 

i — 1 

2.09 

0.49 

0.55 

Ewing  (1957) 

1.11 

2.09 

o.8i 

0.44 

i — 1 

i — 1 

• 

i — I 

2.09 

1.56 

0 . 14 

i — 1 

i — 1 

i — 1 

2.09 

2.64 

0.10 

i 

8 

0.628 

0.51 

1 

8 

1 . 066 

0.17 

1 

8 

1.29 

0 

• 

ro 

0 

1 

8 

1.62 

-0.91 

1 

8 

2.44 

-2.32 

i 

8 

3.76 

OJ 

0 

• 

OJ 

1 

1 

8 

4.94 

1 — 1 

0 

• 

1 — 1 

i 

8 

7.15 

-0.77 

1 

Sezawa  and 

...  ..... 

i 

5 

1.7 

-3.82 

Kanai  (1940) 

1 

5 

3.3 

-1.64 

1 

1 

5 

5.2 

-0.98  j 

1 

20 

2.7 

-3-14 

j 

1 

20 

5.0 

-0.39 

1 

20 

13.2 

-0.99 

1 

00 

2.0 

-2.0 

1 

00 

6.0 

-9 . 6 

1 

00 

9.0 

-14.8 

Positive  Rayleigh  constant  means  retrograde  particle  motion. 
Negative  Rayleigh  constant  means  prograde  particle  motion. 
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on  a  structure  obtained  by  other  means.,  and  the  value  of 
this  check  will  increase  if  it  is  possible  to  determine 
the  Rayleigh  constant  for  a  range  of  different  frequencies. 

The  fact  that  the  Rayleigh  constant  is  sensitive 
to  a  superficial  layer  of  sediments  was  used  by  Oliver  and 
Ewing  (1957)  to  explain  a  discrepancy  between  the  observed 
Rayleigh  constant  and  the  one  predicted  from  a  two  layered 
model . 

A . 3  Necessary  Conditions  for  Sezawa  Waves  to  Exist 

In  section  A.l  it  was  pointed  out  that  for  M2 
waves  a  cut  off  frequency  exists  corresponding  to  a  minimum 
value  of  kH,  i.e.  a  maximum  value  of  L/H.  This  maximum  value 
can  be  expressed  as  a  function  of  \i/\±'  and  vs/vg 1 .  This 
function  was  calculated  by  Yamaguchi  and  Sato  (1955)  and 
their  results  are  shown  on  figure  19.  To  convert  the  wave¬ 
length  L  used  in  the  figure ,  into  frequency ,  a  knowledge 
of  the  phase  velocity  the  waves  is  necessary.  An 

exact  value  of  this  may  be  obtained  from  the  dispersion 
curves  (e.g.  Sezawa  and  Kanai,  1935)  or,  if  applicable,  from 
one  of  the  published  dispersion  curves  for  M2-waves.  An 
approximate  value,  however,  may  be  derived  from  the  following 
information,  given  by  Kanai  (1955)* 

1°  When  the  conditions  for  the  existence  of 
Stoneley  waves  is  satisfied: 

vStoneley  <  Sl2  <  vs 


. 
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MODEL  : 


1 


H 


^  vs  6  Wavelength  L 


T 


x  * p  %  y* 


*  r  vs  & 

(  Poissons  condition  ) 


FIG.  19.  Critical  curvesthat  for  given(-~?)  divides  the  diagram 
into  two  regions  ,  the  upper  one,  where  only  Mj  Rayleigh 
waves  can  exist,  and  the  lower  one,  where  both  Mj  and 


M 2  waves  are  allowed  to  exist. 
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r$ 


#  y  \2 

When(~^“j  K  0.8456  M2  waves  cannot  exist. 


( After  Yamaguchi  and  Sato) 


A  10 


When  the  condition  for  the  existence  of 
Stoneley  waves  Is  not  satisfied: 


v  *  <  CMO  <  v 
s  M2  s 


The  Implication  from  2°,  that  v  *  <  v  ,  is  In 

s  s 

accordance  with  the  fact  that  waves  do  not  exist  for 

vg 1  >  vg  unless  the  conditions  for  the  existence  of  Stoneley 

waves  is  satisfied. 

The  velocity  of  the  Stoneley  wave,  quoted  in  1° 
"falls  between  the  velocity  of  Rayleigh  waves  and  that  of 
transverse  waves  in  the  medium  of  greater  acoustic  velocity" 
(Ewing,  Jardetzky,  and  Press,  1957 .>  p.  112). 


B29820 


